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INTRODUCTION 


The  use  of  physically  small  ^tennas  and  antenna  arrays  can  produce  large  benefits, 
pMicularly  considering  the  limited  space  available  on  missiles,  satellites,  and  aircraft! 
Unfortunately,  if  the  desired  operating  frequency  is  low,  physically  small  also  means 
electncally  small.  For  example,  if  the  operating  frequency  is  around  500  MHz,  the 
wavelength  is  0.6  m,  60  cm,  or  approximately  two  feet.  On  a  four-inch  diameter  missile 
It  IS  clear  Aat  the  largest  antenna  dimension  available  is  liable  to  be  only  a  fraction  of  a 
wavelength,  i.e.  electrically  small.  A  little  reading  of  antenna  books  and  literature 
devoted  to  antenna  theory  and  practice  will  soon  lead  to  the  conclusion  that  electrically 
small  antennas  present  severe  problems  in  terms  of  efficiency,  operating  bandwidth, 
beam-widm,  antenna  patterns,  gain  and  directivity.  Experimental  work  will  quickly 
reinforce  toese  conclusions.  Nonetheless,  the  possible  variations  in  antenna  and  antenna 
a^ay  design  ajppe^  infinite,  so  one  can  always  hope  to  come  up  .with  the  magical 
efficient,  highly  directive,  wide  bandpass  (if  that  is  what  is  desired)  antenna  that 
somehow  “beats  the  game.” 


1  j  need  to  use  electrically  small  antennas  and  the  problems  of  implementation  have 
led  to  a  long  continuing  effort  on  finding  the  theoretical  limits  on  small  antenna 
performance  (References  1  through  11).  That  this  work  has  continued  for  so  long  is 
perhaps,  due  to  the  large  payoff  that  could  be  obtained  with  at  least  acceptable 
electncally  small  antennas,  and  some  doubt  and  confusion  as  to  the  universal  validity  of 
the  precedii^  theoretical  results.  Indeed,  Grimes  (References  12  through  24)  appears  to 
claim  that  Q  and  bandwidth  restrictions,  arrived  at  by  Wheeler  (Reference  2)  Chu 
(Reference  1),  Harrington  (Reference  3)  and  others,  are  much  too  restrictive  and  apply 
only  to  certain  simple  antennas.  He  asserts  that  much  better  performance  can  be  obtained 

by  using  antennas  consisting  of  the  proper  combination  of  electric  and  magnetic 
multipoles.  ® 

It  is  not  the  purpose  of  this  report  to  determine,  once  and  for  all,  what  are  the 
prformance  limts  of  electrically  small  antennas.  It  is  rather  to  lay  a  solid  theoretical 
foundation  and  physical  understanding  of  the  radiation  and  stored  energy  properties  of 
antennas  needed  to  find  these  theoretical  limits.  This  foundation  is  badly  needed 
because,  while  the  papers  quoted  generally  deal  in  some  fairly  involved  mathematics  to 
prove  the  p^icular  points  of  the  particular  papers,  they  are  not  only  difficult  for  the 
reader  to  follow,  but  all  suffer  from  some  lack  of  careful  interpretation  of  the  meanings  of 
Foyntmg  s  theorem,  equivalent  circuits,  and  Q.  For  this  purpose  a  careful  derivation  of 

Poynting’s  theorem  for  both  the  frequency  domain  [exp(j(Dt)]  and  the  time  domain 

[cos(0t  +  <j))]  forms  for  the  electromagnetic  fields  is  initiated  in  the  next  section. 


3 


NAWCWDTP8419 


THE  COMPLEX  POYNTING’S  THEOREM  (TIME  HARMONIC  CASE) 

Technical  books  and  papers  on  electromagnetic  theory,  at  least  in  modem  times, 
generally  assume  electromagnetic  fields  of  the  form 


E  =  Eoexp(y£or) 

(la) 

H  =  Hoexp(7n)r) 

(lb) 

where  Eq  and  Hq  are  complex  vectors. 

However,  such  complex  fields  have  no  physical  existence — it  could  be  said  that 
"Mother  Nature,"  only  knows/contains  real  numbers,  not  complex  or  imaginary.  Often 
electromagnetic  texts  (References  25  through  36)  will  note  that  the  use  of  complex  fields 

is  merely  a  mathematical  convenience  and  that  the  actual  fields  are  of  sin  cot  or  cos  cot 
form,  and  what  is  meant  or  understood  is  actually  the  real  parts  of  Equations  la  and  b. 

Breaking  Eq  and  Hq  into  real  and  imaginary  parts. 


Wl 

o 

il 

Wl 

o 

1 

o 

(2a) 

1 

o 

Iffi 

11 

o 

Iffi 

(2b) 

it  is  simple  to  show  that 


5  =ReE  =  Eo  coscw  +  Eo  sinow 

(3a) 

^  =  Re  H  =  Hq^  coscor  +  Hq  sinew 

(3b) 

It  should  be  noted  that  5  and  ^  are  not  physically  realizable  either,  but  only  in  the  sense 
that  they  have  neither  beginning  nor  ending  in  time.  Otherwise  Equations  3a  and  b  are 
perfectly  general  for  sinusoidal  time  dependence  of  the  fields. 

For  linear  problems  in  electromagnetics,  (i.e.,  problems  in  which  the  field  quantities 

S  and  or  in  circuit  problems,  current  and  voltage,  appear  only  to  the  first  power)  it 
generally  suffices,  and  is  mathematically  simpler,  to  use  the  complex  or  time  harmonic 
form,  work  the  problem,  and,  if  necessary,  take  the  real  part  upon  completion. 
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For  nonlinear  problems,  however,  this  procedure  will  lead  to  erroneous  results.  For 
example,  suppose  there  is  a  device  whose  output  is  proportional  to  voltage  squared,  i.e., 

L  =  KV'^. 


If  the  form  V  =  exp(jQ)t)  is  used,  then 


L  =  KVq  exp{j2cix) 


and  if  the  real  part  is  taken,  L  =  KVq  coslcot. 

On  the  other  hand,  if  the  V  =  Vo  cos  cot  form  is  used. 


=  Vq  cos^  ox  =  Vq  (1 +cos2tBr)/2 


so  that 


L  =  KVi(l  +  cos2<m)/2  , 


which  is  clearly  not  equal  to  the  first  result. 

Much  of  the  work  previously  referenced  is  based,  at  least  in  part,  on  the  Poynting 
vector  and  Poynting's  theorem.  There  are  two  forms  of  Poynting's  vector  in  general  use, 
the  complex  Poynting  vector 


Sc  =  (ExH*)/2 


(4a) 


and  the  time  domain  Poynting  vector 


S  =  giC» 


(4b) 


I^antenna  problems  generally,  H  is  proportional  to  E  (and  g  to  »9^)  so  it  is  clear 
that  Sc  and  S  are  proportional  to  the  fields  squared  and  are,  in  that  sense,  nonlinear.  It 
follows  that  Sc  should  be  used  and  interpreted  with  caution. 

In  most  physics  texts  (Reference  25),  the  time  domain  Poynting  vector  is  interpreted 
as  the  instantaneous  intensity  of  energy  flow  at  a  point  in  an  arbitrary  electromagnetic 
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field,  i.e.,  the  energy  per  second  crossing  a  unit  area  whose  normal  is  oriented  in  the 

direction  of  the  vector,  5  x  »?.  Correspondingly,  the  real  part  of  the  complex  Poynting 
vector  is  interpreted  as  the  mean  intensity  of  the  energy  flow  in  a  harmonic  electro¬ 
magnetic  field. 

It  is  nonetheless  true  that  most  antenna  texts  (References  26  through  36)  choose  to 
work  with  S^,  probably  for  the  sake  of  mathematical  simplicity.  Doing  this  has  led  to 
some  misinterpretation  of  Poynting’s  theorem  and  its  application  to  antennas.  However, 
to  refresh  the  reader’s  memory,  Poynting's  theorem  in  complex  form  is  derived  first. 
(While  this  derivation  can  be  found  in  many  texts,  the  steps  are  important  in 
understanding  what  Poynting's  theorem  does  and  does  not  say.) 

A  useful  vector  identity  is 


V  E;cH*  =  H*V;tE-E  Vj«:H*  .  (5) 

Maxwell's  equations  for  the  time  harmonic  case  are 

V  jc  E  =  -  jcofM  ,  (6a) 

VjcH  =  7  +  7meE  .  (6b) 

Substituting  these  into  Equation  5  gives 

VEA:H*  =  ym(eEE*-/iHH*)-E- J*  .  (7) 


Integrating  both  sides  with  respect  to  volume  gives 


I  V-Exn*dv  =  jO)j  j  E-J*d\  .  (8) 

V  V  V 


Equation  8  is  clearly  a  mathematical  fact,  given  Maxwell’s  equations  and  the  time 
harmonic  form  of  the  fields,  even  though  that  form  is  not  physically  realizable.  It  is 
pertinent,  however,  to  ask,  "Integrate  over  what  volume?"  The  answer,  of  course,  is  any 
volume  one  may  choose.  It  will  usually  be  reasonable  to  choose  a  volume  enclosing 
the  antenna,  although  Equation  8  is  still  true  even  if  this  is  not  done.  Usually  it  simplifies 
the  algebra  to  use  a  coordinate  system  such  that  the  antenna  is  located  at  the  origin 
although  that  is  a  convenient,  not  a  necessary  choice.  Coordinate  systems  and  their 
location  are  the  "inventions"  of  mankind,  and  one  is  free  to  dispose  of  them  as  wished 
although  it  is  usually  preferable  to  choose  them  such  that  can  one  at  least  do  the  algebra. 
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Any  finite  sized  antenna,  no  matter  how  large  physically  must,  by  physical  reasoning, 
look  like  a  point  source  if  the  observer  is  far  enough  away.  Thus  spherical  coordinates 
are  the  natural  coordinates  for  most  antenna  problems,  and  a  logical  volume  in  most  cases 
would  be  a  sphere  with  the  antenna  at  its  center. 

The  next  step  is  to  apply  the  divergence  theorem  to  Equation  8,  ending  up  with 
Poynting's  theorem  for  the  time  harmonic  case,  i.e.. 


(l/2)f  ExH*-da  =  am/2)J  (eE  E*-/tH  H*)dv-(l/2)|  E-/*Jv  .  (9) 

A  V  V 


Equation  9  is  essentially,  a  mathematical  identity,  given  Maxwell's  equations  and  the 
assumption  of  time  harmonic  fields.  However,  it  is  the  "almost"  any  volume  for  the 
volume  integration  on  the  right-hand  side  that  causes  problems. 

The  solutions  to  Maxwell's  equations  in  spherical  coordinates  are  composed  of 
products  of  spherical  Hankel  functions  (actually  of  the  first  and  second  kinds),  their 

derivatives,  associated  Legendre  polynomials,  their  derivatives,  and  functions  of  (])  in  the 

forms  of  exp(jm<|)),  cos  m<|),  and  sin  m(|).  These  are  the  spherical  harmonic  solutions 

(Reference  25). 

Typically  the  antenna  is  considered  as  a  point  source  at  the  origin  (although  it  is  clear 
that  actual  antennas  must  have  finite  extent)  and  often  only  spherical  Hankel  functions  of 
the  second  kind  are  used,  which  are  generally  of  the  form  (Reference  25) 


exp(-7A:r)  X  a„(l/kr)"  , 
n=l 


whereas  the  first  kind  are  of  the  form 

tx^{jkr)f,  an*{\lkr)  , 
n=l 


on  the  basis  that  radiation  propagates  only  outward  from  the  antenna  as  signified  by 
exp(-jkr). 

If  the  antenna  is  not  at  the  origin  or  is  considered  to  be  finite  in  extent,  spherical 
Bessel  functions  would  have  to  be  utilized  as  well,  at  least  inside  some  sphere  containing 
the  antenna  (Reference  35).  Regardless,  the  spherical  Hankel  functions  contain  a 
singularity,  and  become  infinite  at  r=  0. 

The  divergence  theorem  only  applies  to  volumes  that  contain  no  singular  points 
(References  37  and  38).  Thus  the  origin  must  be  excluded  from  Equation  9  if  the 
solution  is  to  contain  spherical  Hankel  functions.  This  can  be  done  by  excluding  a  small 
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sphere  of  radius  ai  around  the  origin  (any  other  non-spherical  volume  around  the  origin 
would  presumably  suffice  but  would  no  doubt  greatly  complicate  the  integration  process). 

Having  excluded  the  origin,  and  hence  the  singular  point  (References  37  and  38),  the 
surface  integral  on  the  left  must  be  performed  over  the  complete  surface,  where  the  origin 
is  surrounded  by  an  arbitrary  inner  sphere  of  radius  a.  and  an  outer  sphere  of  radius  a^ 
which  is  as  large  as  desired.  To  signify  this.  Equation  9  is  rewritten  as 


(l/2)f  BxH*  da 
A 


V  V 


where  the  subscripts  a.  and  ao  indicate  that  the  integration  must  be  over  the  complete 
surface  including  the  inner  sphere  used  to  exclude  the  singular  point  at  the  origin  (see 
Figure  1). 


FIGURE  1.  Complex  Surface  of  Integration  for  the  Complex 
Poynting  Theorem  With  Actual  Antenna  Excluded. 


There  is  no  particular  limit  on  the  size  of  aj.  Frequently  the  current  density  J  at  the 
antenna  is  not  known  and  is  difficult  to  solve  for.  Thus  it  is  common  practice  to  assume  a 
current  density  or  current  on  the  antenna  that  is  "reasonable,"  and  work  from  there.  An 
example  is  the  assumed  sinusoidal  current  distribution  on  a  dipole  antenna.  For  a  half 
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example  is  the  assumed  sinusoidal  current  distribution  on  a  dipole  antenna.  For  a  half 
wave  dipole  using  this  assumed  current  distribution,  "exact"  electromagnetic  fields  can  be 
derived  (Reference  39).  It  is  shown  in  Appendix  A  that  these  "exact  fields"  cannot  be 

correct  in  the  sense  that  tangential  E  does  not  go  to  zero  at  the  antenna  surface,  as  it 
should  for  a  good  conductor.  Thus  there  is  motive  for  making  aj  the  radius  of  the 
smallest  sphere  that  will  just  enclose  the  antenna.  The  volume  between  a.  and  a^  then  has 
J  =  0,  and  Equation  10  becomes 


=  {j(Ol2)\  [eE  E*-tiYi  Yi*)ds  (11) 


in  the  outer  region.  This  is  the  situation  shown  in  Figure  1  and  will  be  the  case  used 
generally  in  the  remainder  of  this  report. 

Aside  from  the  mathematical  convenience  of  not  having  to  deal  with  the  volume 
—  — ^ 

integral  over  E- J  ,  there  is  another,  perhaps  more  cogent  reason  to  exclude  the  actual 
antenna  where  J  exists. 

Consider  Figur^2.  In  this  case  the  integration  includes  the  volume  of  the  actual 
antenna,  V,,  where  J  exists  and  the  fields  are  given  by  Ej  and  Hj .  Exterior  to  V,  in  the 

sphere  of  radius  Uq  is  the  volume  Vj  with  the  fields  %  and  H2.  It  is  H2and  E2  that  can 
be  expresse^in  terms  of  the  exterior  spherical  harmonics  and  thus  have  the  singularity  at 
r  =  0.  Ej,  Hj ,  and  J  need  have  no  such  singularity  and  physical  argument  can  be  made 
that  since  Vj  contains  the  actual  antenna,  no  such  singularity  exists. 

Equation  8  can  now  be  written  as 


(1/2)/  Ex'R*  dd 

A 


j  VEi  xHi*dv+  J  V-E2xH2*dv  =  jwJ  (ciEi -Ej  Hi  *Wv 

v;  y, 


+ JG>  /  (eoE2  •  ^2  *  -^0^2  •  H2  *)</v  -  /  EyJ*dv 

y,  y, 


which  is  true  without  any  particular  qualification. 
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FIGURE  2.  Complex  Surface  of  Integration  for  the 
Complex  Poynting  Theorem  With  Actual  Antenna 
Included. 


On  the  other  hand,  if  one  tries  to  work  with  Equation  9,  noting  the  lack  of  singularity 
so  that  an  inner  and  an  outer  sphere  need  not  be  considered,  one  would  have 


f  E2xH2*  da 
A 


—  j(0  !  (^lEi  •  Ej  *  •  Hj  *^d\ 

V 


+  ja)j  (eoE2  E2*-AtoH2  H2*yv- J  Ei  J*t/v 


Now,  while  E2,  Hj,  E2,  and  H2  are  all  solutions  to  Maxwell's  equations  and  are 
constrained  by  the  usual  boundary  conditions,  they  are  not  the  same  functions,  in  general, 
as  shown  by  the  lack  of  singularities  in  Ej  and  Hj  and  the  existence  of  J  in  only  V,. 

However,  if  E  and  H  (or  at  least  V  ExH*)  are  not  the  same  functions  of  the 
coordinates  in  both  V,  and  V^,  the  divergence  theorem  does  not  apply. 

Thus  if  we  are  considering  a  volume  with  /(r)  in  one  part  and  g(r)  in  another  part, 
the  divergence  theorem  can  be  applied  to  the  surface  of  that  part  of  the  volume  that 
contains  /(r)and  to  the  surface  of  that  part  which  contains  |'(r)  separately,  but  not  to  the 
surface  that  encloses  both  volumes  together  (References  37  and  38).  The  volume  must  be 
unifunctioned. 
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The  complex  Poynting  theorem  is  usually  written  in  the  form  of  Equation  9  and  it  will 
be  observed  that  the  surface  integral  in  Equation  9  is  just  the  power  flow  across  a  surface 
surrounding  the  antenna  and  the  real  part  of  that  integral  is  the  time  average  or  radiated 
power  flow  across  that  surface.  This  is  true  if  the  surface  integral  is  only  evaluated  at  a„ 
but  this  has  little  to  do  with  Poynting's  theorem  per  se;  i.e.,  this  could  have  been  simply 
stated  without  ever  deriving  the  theorem. 

In  free  space  E  and  p  are  real,  thus  the  first  volume  integral  on  the  right  hand  side  of 
Equation  9  is  imaginary.  Since  aj  is  taken  to  be  the  radius  of  a  sphere  large  enough  to 

enclose  the  antenna  (where  J  is  zero  outside  such  a  sphere),  the  second  volume  integral 
on  the  right  is  zero,  leaving  the  right-hand  side  of  Equation  9  imaginary.  It  follows  that 
the  left-hand  side  of  Equation  9  must  also  be  imaginary,  as  indeed  it  is  if  the  integration  is 
over  the  complete  surface  as  in  Equation  10. 

Physical  reasoning  shows  that  this  must  be  the  case.  If  the  surface  integral  over 
(l/2)Re(ExH*)  is  the  time  average  radiated  power,  and  then  it  must  be  the  same  for 

any  sphere  surrounding  the  antenna,  i.e..  Ref  E  x  H  *-da  is  independent  of  the  spherical 


However,  performing  the  integral  over  the  complete  surface  is  the  same  as 


f  Exli*  dd 

=  f  Exli*Nida 

+  f  ExK* 

.A 

a,  A 

But  the  outward  normal  at  a^  is  Ng  =  r,  and  at  aj,  the  inner  normal  is  iV,-  =  -f  ,  so  the 
real  parts  must  cancel. 

As  a  corollary,  since  Re[E  x  H  *  -da]  is  independent  of  radius  and  the  increment  of 
area  is 


da  =  r^sin0  dd  d(j) 


(12) 


the  real  part  of  E  x  H  *  must  go  as  1/r^  so  the  integration  will  be  independent  of  r,  and 
only  those  portions  of  E  and  H  that  go  as  1/r  contribute  to  the  radiated  power. 

To  sum  up  this  section,  those  who  write  Poynting's  theorem  (Reference  26)  in  the 
form  of  Equation  9,  then  state  that  the  real  part  of  the  left  hand  side  of  Equation  9 
represents  radiated  power  are,  at  best,  misleading.  They  should  state  that  the  time 

average  radiated  power  is  indeed  given  by  (l/2)Ref  E  jc  H  *  da  on  a  sphere  external  to 

A 

the  antenna,  but  that  this  is  not  actually  part  of  Poynting's  theorem  which  is  correctly 
expressed  by  Equation  10  (with  care)  or  by  Equation  1 1  for  the  case  of  antennas  where 
there  is  a  singularity  at  the  origin. 


11 


NAWCWDTP8419 


Finally,  Poynting's  theorem  can  be  confusing  in  another  way,  insofar  as  the  dielectric 

constant,  e,  and  magnetic  permeability,  p,  are  concerned.  At  least  one  text  (Reference 
30)  notes  that  the  imaginary  part  of  the  first  integral  on  the  right  hand  side  of  Equation  9 
(which  should  be  written  as  Equation  10),  upon  being  multiplied  by  j,  represents 

dissipated  energy.  If  one  understands  this  to  mean  that  e  is  the  complex  dielectric 
constant  as  usually  def  ned  from 


VxH  =  Ic  +  BDIdt  =  oE  +  jmeofi'E  =  j(oeQ[e' - jcT/toeo)^  =  , 


one  would  be  wrong.  In  deriving  Poynting's  theorem  the  conduction  current  was  already 

taken  into  account.  Thus,  in  Equations  8, 9,  and  10,  e  is  real  (=  z^'),  unless  some  kind  of 
losses  other  than  conduction  losses  have  been  included  in  the  displacement  current  term 
in  Maxwell's  equation. 

That  this  is  so  can  be  seen  as  follows.  Consider  some  volume  in  space  characterized 
by  8,  p,  and  a  where  there  is  a  current  density  J,  and  electric  and  magnetic  fields  E  and 

H.  It  is  assumed  that  in  this  volume  these .  quantities  are  everywhere  finite  (no 
singularities“it  can  be  argued  that  for  any  physically  real  situation  this  must  be  the  case). 
In  tiiis  case.  Equation  9  applies  without  any  need  for  an  inner  and  an  outer  surface — only 
the  outer  surface  need  be  used. 

It  follows,  in  this  case,  that  the  left-hand  side  of  Equation  9  will  have  both  real  and 
imaginary  parts,  with  the  real  part  representing  the  radiation  out  of  the  volume  V.  It  will 

be  assumed,  for  this  example,  that  p  =  Po ,  and  the  current,  J,  is  given  by  Ohm's  law. 


J  =  oE 


(13) 


where  a  is  real. 

If  it  were  assumed  that  is  the  complex  dielectric  constant  as  given  by 
Eg  =  80(8'  -ja/co  Eq)  with  e'  and  <5  real,  the  right-hand  side  of  Equation  9  becomes 


/mj  EQ(e'-jal(oeQ)E-B*dv-jo)\  poH  H*Jv-J  oE  E*Jv 
V  V  V 


=  7©/  [eos'E  E*-fioH-li*)dw 
V 


That  is,  the  right-hand  side  of  Equation  9  would  be  totally  imaginary,  which,  of 
course,  is  impossible  since  power  is  radiated  and  the  left-hand  side  has  a  real  part. 
Clearly  one  has  to  be  very  careful  in  interpreting  and  applying  the  complex  Poynting's 
theorem. 
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TIME  DEPENDENT  POYNTING'S  THEOREM  (GENERAL  CASE) 


Beginning  with  the  vector  identity 

V  •  S{t)  X  ^(t)  =  ^(t)  •  VxF(t)  -  S(t)  •  Vx  ^(f)  (14) 

where  S(t)  and  f?^(t)  are  general  real  (physical)  time  dependent  fields, 

Maxwell’s  equations  are  now 

V xS{t)  =  -dS{t)ldt  ,  (15a) 

Vx^{t)  =  ^{t)  +  ^{t)ldt  .  (15b) 

Substituting  Equation  15a  and  b  into  Equation  14  gives 


Integrating  both  sides  over  a  volume  in  accordance  with  the  previous  section  and 
applying  the  divergence  theorem: 

^  S{i)x^{t)-dd  =  I  V •  S{t)x1i^{t)d\ 

A  V 

(17) 

=  d8{t)l dt  +  S{t)- dD{t)l dt\dw -  \  S{t)  f{t)dw  . 

V  y 

Assuming  that  f  (t),  »?^(t),  and  /  (t)  has  the  sinusoidal  time  dependence  of  Equation 
3a  and  b  and  that 


V  =eS 

(18a) 

€  -  , 

(18b) 

Equation  17  becomes 
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^  Sx1ii-da=-\  [}j5¥  d»ldt  +  eSdSldi)d\-\  g  pdv  .  (19) 

A  V  V 


Equation  (19)  can  just  as  well  be  written  in  the  form 


The  same  rales  regarding  the  divergence  theorem  apply  as  in  the  previous  section  that 
the  volume  must  contain  no  singularities  and  that  S,  9-,  and  ^  must  be  the  same 
functions  throughout  the  volume. 

Choosing  an  inner  sphere  that  encloses  the  antenna  and  any  outer  sphere  to  meet  these 
criteria, 


.A 


\  [\}e»-1¥l2-\-eS-Sll)dv  . 
V 


(21) 


One  can  recognize  the  sum  of  e(5  •  5)/2  and  p(  as  the  total  electric  and 

magnetic  energy  per  unit  volume  in  V.  Thus 


f  Sx'^‘dE 


=  -dUldt 


(22) 


where  U  is  the  total  energy  in  the  volume,  V. 

What  is  the  relationship  between  Equations  21  and  11?  After  all,  complex  numbers 
cannot  truly  represent  physical  fields.  From  Equation  3a  and  b 


g‘S  =(Eo  coscw  +  Eq  sintat|  |Eo  costat  +  Eo  sintorj 
=  (jEo)^cos^dar  +  (£o)^sin^tar+2Eo  Eq  sin  tacos  tor 


Finally,  with  the  use  of  trigonometric  identities. 
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5  •  5/2  =  (l/4)((£6f  +(^6f  )  +  (l/4)((56f  -{E^f)cos2QX 
+  (1  /  2)£'q  •  Eq  sin  2cot 

and  it  follows  that 

^  •  ^/2  =  (l/4)((H6f  +  (H6f )  +  (l/4)((H6f  -  )cos2g» 
+  {1/2)Hq  -Hq  sm2ciX  . 

Now 

f  •  I*  =  (^'  -  =  (E(,f  +  (E^f 

and 

=  +(^of  • 


(23a) 


(23b) 


(24a) 


(24b) 


Therefore,  e  5  •  5  *  /  4  is  the  time  average  energy  per  unit  volume  in  the  electric  field, 
and  \lH  -H  */4  is  the  time  average  magnetic  energy  per  unit  volume. 

On  the  other  hand 

-  ^£9  costat  +  Eq  sin  costor  +  Hq  sin  tar j 

=  (l/2)|£'o  ^^0  ■’■^0  ^^0  )  +  (1/2)|£o  ^Hq  -Eq  xHq  jcos2taf  (25) 
+  (l/2)|£'o  ^Hq  +Eq  xHq  jsin2<at 

while 
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ExH*  =  [EQ-jE^'y^x(Eo+jEQy 

Eq  xHq  +£0  j  +  j|£o  ^^0  ~^0  ^^0  j 


(26) 


Comparing  Equations  25  and  26  one  finds 


{l/2)Re(ExH*)  =  (Sx‘») 


(27) 


where  the  bar  means  time  average.  The  imaginary  part  of  E  x  H*,  however,  bears  no 
obvious  relation  to  any  part  of  Equation  25  except  that  it  is  the  difference  instead  of  the 

sum  of  the  quantities  in  brackets  in  the  coefficient  of  sin  2(0t  in  Equation  25. 


SUMMARY  AND  DISCUSSION  OF  POYNTING’S  THEOREMS 


Given  time  domain  variation  of  the  fields  and  currents  (Equation  3a  and  b),  the  time 
domain  form  of  Poynting's  theorem  is  , 


j  Sx‘»-da=~!  (ii'»-^/2  +  £S-S/2)dv-  j  S  fdv 
A  V 

where  the  volume  integrated  over  can  be  any  volume  with  the  surface  integral  being  over 
the  complete  surface  of  that  volume.  This  is  provided  the  volume  integrated  over 

contains  no  singularities  and  S,  and  /  are  the  same  vector  functions  throughout  that 
volume. 

The  time  harmonic  or  complex  form  of  Poynting's  theorem  is 


(1/2)^  ExH*da=j2(o\  [eE  E* lA-pJti  H* IA)dv-{\l2)\  E  Jdv 
A  V  V 


This  is  actually  Equation  9  rewritten  slightly  to  emphasize  the  time  average  relations 
between  Equation  9  and  the  previous  equation.  The  same  conditions  as  to  the  volume, 
surface,  and  functions  integrated  over  apply,  and  it  should  be  noted  that  while  Equation  9 

is  mathematically  correct,  E  and  H  are  not  physically  realizable  fields. 

The  usual  solutions  (the  spherical  Hankel  functions)  for  radiation  in  spherical 
coordinates  contains  a  singularity  at  the  origin,  therefore,  if  these  fields  are  to  be  used  the 
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volume  must  exclude  the  origin.  Real  antennas  are  finite  in  extent  and  generally  the 
fields  and  currents  in  the  volume  actually  containing  the  antenna  are  not  functionally  the 
same  as  the  fields  and  currents  (usually  zero)  exterior  to  the  antenna,  even  though  these 

fields  must  meet  the  usual  boundary  conditions  of  continuous  tangential  S  and  ^  with 
the  interior  fields.  Therefore,  the  volume  must  exclude  the  actual  antenna  as  well  as  the 
origin  for  Poynting's  theorem  to  be  applicable. 

Typically  /  is  zero  exterior  to  the  antenna  and,  given  that  we  are  usually  interested  in 
a  volume  that  surrounds  the  antenna  in  order  to  determine  radiated  power,  Poynting's 
theorem  reduces  to 


f  Sxl^-da 


A 


[\)S¥-‘»l2-\-eS'Sli)dv 


(28) 


where  the  volume  V  is  the  volume  between  any  inner  surface  Aj,  which  surrounds  the 
antenna  and  any  outer  surface,  A^,,  which  also  surrounds  the  antenna. 

The  equivalent  form  for  the  complex  Poynting  theorem  is 


2 


_A 


=  ;2mJ  (eE  E*/4-|iH-H*/4Wv 
V 


(29) 


Both  sides  of  Equation  29  are  purely  imaginary. 

Because  closed  form  integration  over  anything  but  spherical  surfaces  in  spherical 
coordinates  can  be  at  least  tedious  if  not  impossible,  Aj  and  A^  are  easiest  to  work  with  if 
they  are  spherical  so  that  a;  is  the  radius  of  the  smallest  sphere,  which  will  contain  the 
antenna,  whether  the  antenna  is  spherical  or  not,  and  a^  is  the  radius  of  any  sphere 
exterior  to  this  with  the  same  origin,  in  which  case  Poynting's  theorems  become: 


f  Sx^  -da 


lA 


I  ’  0 


J  (\L^  n  +  es  ■  S  I2)dw  =  -du ! dt 

dt  Y 


and 


ExH*dd 


=  j2(o\  [eE  E*IA-\^i-H*IA)dw 
V 


Note  that  eE-E* lA  is  the  time  average  of  eS  'S 12,  i.e.,  the  time  average  of  the 
total  electric  energy  per  unit  volume  in  V.  The  same  relation  holds  for  \iH  H  *  !  A  and 
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The  right-hand  side  of  Equation  21  is  the  time  derivative  of  the  total 
electromagnetic  energy  in  the  volume  V,  whereas  the  right-hand  side  of  Equation  1 1  is 

j2©  times  the  difference  between  the  time  averages  of  the  magnetic  and  electric  energies 
inV. 

Finally,  it  is  noted  (and  this  is  separate  from  Poynting's  theorem)  that 


Pa  =  (l/2)Re^  ExH*  dd  (30) 

A 


where  the  above  integral  is  around  any  surface  surrounding  the  antenna  and  P^  is  the  time 
average  power  flow  across  that  surface,  i.e.,  the  radiated  power,  while 


^  Sx^  da 
A 


(31) 


is  the  instantaneous  power  flow  across  that  surface. 
Note  that,  in  general 


J  [fj:»-‘»/2  +  eS-S/2)d\ 


V 


and 


j2G)j  (eE  E*/4-\iH  H*/4)dw 
V 


as  the  right-hand  side  of  the  latter  doesn't  usually  exist  in  most  cases  of  interest. 

A  reasonably  clear  picture  of  Poynting's  theorem(s),  power  flow  across  a  surface 
surrounding  an  antenna,  and  energy  contained  in  the  volume  surrounding  the  antenna  is 
now  achieved. 

The  question  arises  as  to  how  to  employ  these  concepts  to  tell  something  about  the 
antenna  Aat  one  needs  or  wants  to  know.  Of  course,  to  carry  out  the  integration,  dot  and 
cross  products  etc.,  the  presumption  is  that  one  already  knows  the  mathematical  forms  of 
the  electric  and  magnetic  fields.  If  this  is  so,  one  already  has  the  information  needed  for 
the  antenna  pattern,  gain,  beamwidth,  sidelobe  levels,  polarization,  etc.  In  addition  to 
these  things  it  would  generally  be  useful  to  know  some  other  properties,  like  input 
impedance  and  bandwidth,  i.e.,  the  behavior  of  the  input  impedance  with  frequency. 
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Engineers  typically  relate  input  impedance  and  bandwidth  to  concepts  such  as 
resistive  and  reactive  power  flow,  equivalent  circuits,  and  the  quality  factor,  Q.  The  latter 
is  usually  found  through  ad  hoc  relations  such  as 


e=/o/Af 


(32) 


where  fo  is  the  operating  or  resonant  frequency  and  Af  the  3  dB  bandwidth  around  f,,.  The 
equation  for  the  Q  is  often  given  also  as  a  relation  between  stored  energy  and  dissipated 
power  as 


Q  =  (oUsIWl 


(33) 


where  is  the  time  average  power  loss  and  is  the  energy  stored  in  the  circuit.  The 
implication,  of  course,  is  that  "high  Q"  circuits  have  a  narrow  operating  bandwidth.  The 
power  lost,  Wl,  can  be  identified  as  the  power  radiated.  Equation  30,  in  the  absence  of 
additional  ohmic  losses  in  the  antenna.  The  reactive  energy  stored  can  be  found  by  the 

proper  integration  over  S  •  B  and  using  the  reactive  fields. 

Papers  on  the  fundamental  limits  of  antennas  have  followed  this  general  procedure  in 
one  fashion  or  another  to  show  that  for  electrically  small  antennas  the  ratio  of  the  reactive 
energy  stored  in  the  antenna  fields  to  the  average  power  radiated  is  large  (References  1, 
3,  and  7).  It  is  then  concluded  with  reference  to  Equation  32  that  since  the  Q  is  high, 
these  are  intrinsically  narrowband  devices.  Although  there  is  obviously  a  certain  validity 
to  this  general  approach,  one  has  to  be  extremely  careful  not  to  push  it  too  far.  Although 
Equations  32  and  33  may  be  wonderful  generalizations  which  seem  almost  intuitive 
through  long  acquaintance,  they  need  careful  examination  when  applied  to  a  particular 
antenna  or  antenna  system.  As  ^ways,  the  devil  is  in  the  details. 

It  is  common  practice  in  antenna  texts  to  use  the  term  "radiation  resistance,"  usually 
found  by  employing  Equation  30  together  with 


PA=llRrl2  (34) 


where  Iq  is  the  input  current  to  the  antenna  and  R  is  the  radiation  resistance.  When  this 
definition  is  applied  to  an  infinitesimal  dipole  one  gets  the  usual  formula  (References  26 
through  28) 


Rr  =  mK^{^zlXf 


(35a) 


or,  for  a  short  dipole  with  a  physically  reasonable  current  distribution 
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Rr  =  2QK^(hlXf- 


(35b) 


where  Az  and  h  are  the  dipole  lengths. 

While  the  radiation  resistance  is  a  seemingly  straightforward  concept,  it  too  can  have 
its  problems.  It  is  one  thing  to  apply  it  to  an  antenna  with  a  single  pair  of  input  terminals 
where  the  current  is  known  at  those  terminals;  it  is  another  thing  to  find  the  meaning  in 
the  case  of  an  antenna  array  (Reference  40). 


RESISTIVE  AND  REACTIVE  POWER  FLOW  AND  EQUIVALENT  CIRCUITS 
Consider  the  series  RLC  circuit  shown  in  Figure  3  with  an  applied  voltage 

V=Vocosu»  .  (36) 


FIGURE  3.  Series  RLC  Circuit 
With  Applied  Voltage,  V. 


The  current  through  the  circuit  must  therefore  be 

/  =  Re(yoe-'‘“/2)  (37a) 

where  the  impedance,  Z,  is 

Z  =  R  +  jX  =  R  +  j{(oL-\l(oC)  (37b) 
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or 

Z  =  |Z|exp(jtan"Hx//?))  =  |Z|e-'^  . 
The  current  then  becomes 


I  =  {Vo  /lZ|)Re  =  (Vo  /|Z|)cos(n»-  0)  =  Iq  cos(<af  -  <j>) 


and  the  power  into  the  circuit  is  thus,  using  Equations  36  and  38, 
P  =  rv  =  /|Z|jcosn»cos(n»  -  0)  . 

Using  the  double  angle  formula, 

P  -  |Vo  /|Z|jcos(af(cos^cos£or  +  sin^sintof)  . 


However, 


(37c) 


(38) 


(39) 


(40) 


cos<l>  =  IZ-^l  +  tan^^  =  costan  ^  —  =  R/|Z| 

R 


and 


(41a) 


sin^  =  tan0/'^l  +  tan^0  =X/|Z|  . 
Substituting  these  into  Equation  40  we  have 


(41b) 


P  =  (Vo^ /2|Zp )[/?(!  +  cos2mf)  +  Xsin2cof]  . 


(42) 


The  first  term  of  Equation  42  ranges  between  zero  and  VqRI\Z'^  =  IqR  and  would 
seem  logically  to  be  associated  with  the  time  dependence  of  the  resistive  power  loss, 

while  the  second  term  ranges  between  -loXH  and  IqXH  and  could  be  associated  with 
the  reactive  power  being  alternately  supplied  by  and  returned  to  the  source. 
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Such  an  interpretation  sounds  plausible  (Reference  24),  but  is,  in  this  case,  technically 
incorrect.  This  is  readily  shown  by  using  the  reactive  power  as  the  derivative  with 
respect  to  time  of  the  reactive  energy,  i.e.. 


Pxi  =  dUxildt 


(43) 


The  reactive  energy  in  the  inductance  is 


Vii  =  f'LI2 


(44a) 


and  in  the  capacitor 


Ua^cvin 


(44b) 


where  Vc  is  the  voltage  across  the  capacitor. 
Now 


Uu  =  (Vo^  /2|Zp)Lcos^(G}f -  ^)  =  (Vo^L/4lZp)  [l  +  cos2(a» -  ^)]  .  (45) 

The  voltage  across  the  capacitor  is 

Vc  =  Re{I/jcoC)  =  Re(-7yo  /  (oC\Z\)txp[j{(ot  -  (p)] 

or 

Vc  =  {Vo/(aC\Z\)sin{(m-^)  .  (46a) 

Thus 

Ua  =  (^0^  /2n}2qzp)sin^(ta  -(j>)  =  /4u)^qzp)[l  -  cos2(g»  -  (p)]  •  (46b) 

The  total  reactive  energy  is  then 
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Uxi  =  [(©L+l/Q)C)  +  (©L-l/foC)cos2(G»-0)] 


(47) 


which  has  an  interesting  symmetry. 

Substituting  Equation  47  into  Equation  43 

Pxi  =  [-vixi2\Z\^)s,m2{ax-(t>)  (48) 


which  is  of  a  different  phase  and  sign  than  {VqX ! 2\2!^^sm2cot  as  could  be  gathered 
from  Equation  42  although  the  magnitudes  of  Equation  42  and  Equation  48  are  the  same. 

One  can  obtain  the  correct  answer  by  finding  the  resistive  power  loss  directly  as 


=  i^R  =  (Vo^iJ/|Zp)cos^(mt-  =  (Fp/?/2|Zp)[l +cos2(a)r-  0)]  (49) 


and  postulating  (note  that  we  cannot  set  =  f X)  that  the  total  power  flow  must  be 
Equation  49  combined  with  Equation  48: 


P  =  (7o^/2|Zp){/?[l  +  cos2(mt-0)]-Xsin2(mr-0)}  .  (50) 


It  only  remains  to  show  that  Equations  50  and  42  are  equal.  They  are  and  this  will  be 
shown  in  Appendix  B. 

The  reactive  power  flow  could  have  been  derived  directly  by  observing  that 


Px  =  P^X  (51) 


where 


Vx  =  Re[i/X]  =  Rt[jXVoej^/z\ 
or 


(52) 


Px  =  -(yo-S^Z|)sin(a)f  -  <f>)(Vo  /|Z|)cos(mf  -  0)  =  -(Vo^X/2lZp)  sin2(mf-0) 


which  is,  of  course,  Equation  48. 
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In  general,  as  shown  in  Appendix  B,  if  the  applied  voltage  and  current  are  given  by 

V{t)  =  Vq  cos(£0f  -  a)  (53a) 

and 

I{t)  =  Iq  cos{cot  -  p)  (53b) 

then  P  =  rV  can  be  written  in  two  ways,  either 

P  =  (Vq/o  /2){cos(a  -  )5)[l  +  cos2(60f  -  j8)]  +  sin(a  -  P)sm  2{cot  -  j3)}  (54a) 


or 


P  =  {VQlo/2)[cos{a- P)[l  +  cos2{cDt-a)]-sm{a- P)sin2{cix-a)]  .  (54b) 


In  any  given  case,  only  one  of  these  two  expressions  can  be  (directly  and  correctly) 
broken  into  resistive  and  reactive  power  flow.  The  importance  of  this  will  become 
clearer  when  resistive  and  reactive  power  flow  in  relation  to  the  Poynting  vector  is 
examined. 

Before  this  is  done,  however,  there  are  a  couple  of  points  to  be  made.  If  one 
considers  a  single  antenna  as  a  (one  port  or  two  terminal)  black  box  which  contains  any 
series/parallel  combination  of  capacitors,  inductors,  and  resistors,  and  if  one  applies  a 
signal  at  any  single  frequency,  one  can  measure  an  input  impedance  consisting  of  a  real 
resistance,  and  a  series  reactance.  The  exact  relation  to  the  various  resistors  and 
reactance  (capacitors  and  inductors)  actually  in  the  box  depends  on  their  values  and 
series/parallel  arrangement.  In  fact,  from  a  measurement  at  a  single  frequency  with  no 
other  a  priori  knowledge  as  to  the  contents  of  the  box,  one  has  no  information  as  to  what 
the  box  actually  contains  other  than,  if  there  is  a  real  part  to  the  input  impedance,  there 
must  be  at  least  one  resistor  in  the  box,  and  if  there  is  an  imaginary  part,  at  least  one 
reactive  element.  Even  if  one  has  a  priori  knowledge  of  what  the  box  contains,  for 
example,  a  series  RLC  circuit,  a  measurement  at  a  single  frequency  suffices  only  to  give 
the  vdue  of  the  resistance  and  the  reactance,  but  not  L  or  C  individually.  Clearly,  if  the 
box  contains  frequency  independent  resistors,  capacitors  and  inductors,  and  one  wishes  to 
find  their  values  and  their  series/parallel  arrangements,  one  must  measure  the  input 
admittance  or  impedance  as  a  function  of  frequency  and  endeavor  to  find  that 
combination  of  frequency  independent  elements  that  gives  the  measured  frequency 
dependence. 

The  alternative  would  be  to  assign  frequency  dependent  elements  such  that  the 

observed  frequency  dependence  is  Z(©)  =  R((a)  +  jX((D).  Regardless,  the  power  flow  at 
any  frequency  can  be  written  in  the  form  of  either  Equation  42  or  Equation  50. 

To  illustrate  this  point,  suppose  the  black  box  contains  not  the  series  RLC  circuit  of 
Figure  3,  but  the  parallel  RLC  circuit  of  Figure  4.  The  input  admittance  is 
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Y  =  G-^jB  =  \IR+j{(oC-\l(oL)  .  (56) 


FIGURE  4.  Parallel  RLC  Circuit  With  Applied  Voltage,  V. 

The  input  impedance  is 

Z  =  \IY  =  {G-jB)l{G^+B^^  .  (57) 

The  power  dissipated  in  the  resistance  is 

Pf^  =  I^R  =  V^  IR  =  V^Gcos'^  cot  =  {v^G  12^1  + coslcot)  .  (58) 

The  capacitor  current  is 

4  =  Re[Voej^/Xc]  =  Re[jcoCVoe^^'\  =  -VoCoCsmciX  (59) 

and  for  the  inductor  current 

1 1  =  RelVo^"^^  !  y^iLj  =  (Vq  /  sin  cot  .  (60) 

The  power  flow  to  the  capacitor  is 

Pc  =  IcV  =  -{v^aC/iymloot  (61) 

and  to  the  inductor. 

Pi  =  IiV  =  {Vq  /2n)Ljsin2a»  .  (62) 
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The  total  power  flow  is  then 


P  =  Pr+Pc  +  Pl  =  {v^  l2)[G{\  +  coslcot)-{a)C-\l(£L)sm2Gx]  (63a) 

or 

P  =  (v{?/2  )[G(1  +  COS  lojt)  -  5  sin  2£»]  .  (63b) 

In  terms  of  the  previous  discussion,  Equation  63  must  be  the  same  as  Equations  42  or 
50,  which  are  identities.  Choosing  Equation  42  as  the  simpler  case,  we  must  have 


{V^I2^G{\  +  cos2mt) - Bsinlm]  =  /2|Zp)[P(l  +  cos2daf)  +  Xsin2taf]  . 


Substituting  from  Equation  57  this  is  the  same  as 


G(1  +  cos2mf)  -  Bsin2<»f  =  [G(l  +  cos2<»t)  -  Bsin2mf]^|Z|^(G^  + 


Since  Z=  IfY, 


|Z|^  =  1/|7|2  =  i/(g2+52J  ^ 


(64) 


one  arrives  at  G(l+  cos  2cot)  -  B  sin  2cot  =  G(l+cos  2m)  -  B  sin  2cot,  which  is  the 
expected  result. 

If  the  black  box  contains  an  antenna,  then  an  antenna  is  by  no  means  an  RLC  circuit 
but  rather  a  distributed  one.  Nonetheless,  at  the  input  terminals  one  will  measure  some 
input  impedance  with  a  power  flow  which  can  be  written  in  the  form  of  Equation  42  or  its 
identity.  Equation  50.  Knowing  the  input  voltage  and  current  one  will  know  whether  the 
input  is  capacitive  or  inductive  and  whether  Equation  42  or  50  is  appropriate  from  the 
aspect  of  reactive  and  resistive  power  flow.  If  one  measures  the  input  impedance  as  a 
function  of  frequency,  one  can  attempt  to  establish  an  equivalent  circuit  for  the  antenna, 
either  a  circuit  that  involves  frequency  dependent  elements  or  a  circuit  involving 
frequency  independent  elements. 

It  is,  however,  a  goal  of  antenna  theorists  to  establish  the  input  impedance 
theoretically.  logically,  the  first  step  is  to  find  the  "radiation  resistance"  via  the  power 
radiated,  i.e.,  the  surface  integral  of  the  Poynting  vector  over  a  sphere  surrounding  the 
antenna.  Here  the  real  part  of  the  complex  Poynting  vector  can  be  used  or,  if  the 
(physical)  Poynting  vector  is  employed,  those  terms  with  radial  dependence  of  l/i^. 
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The  power  radiated  is  generally  found  in  terms  of  the  assumed  antenna  current 
squared,  and  hence  is  related  to  the  term  I^Rj,  where  I  is  the  input  current  to  the  antenna 

and  Rj.  is  the  radiation  resistance.  Employed  this  way,  is  obviously  a  series  element, 
and  will  generally  be  found  to  be  frequency  dependent.  For  small  dipoles  (References  26 
and  28)  die  radiation  resistance  is  given  as 

Rr  =  %07p'{hlXf  (65a) 


for  an  assumed  uniform  current  distribution  and 
Rj.  =  20n^{hlXf' 


(65b) 


for  an  assumed  sinusoidal  current  distribution. 

In  either  case,  R^  goes  as  one  over  the  wavelength  squared,  i.e.,  as  the  frequency 
squared.  Thus  any  series  equivalent  circuit  employing  the  radiation  resistance  involves 
frequency  dependent  circuit  elements.  Thus  one  can  see  the  motivation  in  using  this 
same  integral  to  find  the  radiation  reactance,  and  perhaps  extending  it  to  the  concepts  of 
Q  and  bandwidth. 

Finally  one  can  note  that  while  the  concepts  of  radiation  resistance,  radiation 
reactance,  and  antenna  input  impedance  seem  conceptually  clear  for  antennas  with  a 
single  input  port  (actually  there  are  caveats  even  here),  these  concepts  become  much 
more  complicated  in  dealing  with  antenna  arrays  (antennas  with  multiple  input  ports). 
This  has  been  covered  to  a  limited  extent  in  an  earlier  report  (Reference  40). 


POWER  FLOW  AND  POYNTING'S  VECTOR. 


One  can  return  to  Equation  25  and  rewrite  it  as  (References  14  and  15) 


S  =  S  =  Nq+Ni  coslox  +  N2  sin2G)f 


where 


, 

xHq  -Eo" xHo"]/2  , 


(66a) 


(66b) 

(66c) 


and 


NAWCWDTP8419 


^2  =  ^0  ^^0  ■•'^0  ^^0  • 


(66d) 


If  one  integrates  5 ,  as  in  Equation  31,  on  a  surface  surrounding  the  antenna,  the  result 
will  be  the  instantaneous  power  flow  across  that  surface.  If  the  surface  chosen  is 
spherical  and  centered  on  the  antenna,  then  da  is  given  by  Equation  12  and  da  is  in  the 
direction  of  f.  Thus  only  the  radial  component  of  Equation  66  enters  into  the  integration 
and 


P  =  ^  N^j-r^sindddd^  +  coslm^  N\j.r^sxndddd(p-\-sm2ox^  N2rr^ sin 6ddd<t>  (67) 
AAA 


which  can  be  written  as 


P  =  Mq+Mi  coslox  +  M2  sin  lax  (68a) 

where 

Mq  =  I  NQ^r^  sin6d6d(t>  (68b) 

A 

Mi  =  j  Nif.r^ sinddOd^  (68c) 

A 

M2  =  j  N2rr^ sm6ddd(j)  .  (68d) 

A 

Generally  speaking,  the  complex  time  harmonic  solutions  for  the  fields  in  antenna 
problems  have  the  form  of  F(r,  6,  <|))  exp(-jkr)  exp(j©t) .  It  is  thus  natural  to  combine  cot 
and  -kr  as  cot  -  kr,  which  can  be  written  as  cot^  where 

tr  =  t-{k/(o)r  =  t-r/c  .  (69) 


t,  is  called  the  retarded  time  and  it  simply  means  that  physically  the  propagation  time  may 
be  small  but  it  is  not  zero.  Some  authors,  notably  Grimes,  generally  specify  the  use  of 
retarded  time  in  expressing  power  flow  (Reference  22). 

In  fact  one  can  always  write  Equations  la  and  lb  in  the  form 
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E  =  Eie.x^(jmr)  =  {Ei  -jEi  jexp(7G»r)  . 

(70a) 

H  =  HiexT^jaXr)  =  {Hi  -jHi  jexp(;£Wr) 

(70b) 

so  that 

^  _ /  _ // 

S  =  Ei  cosaXf  +  Ei  sinaXr 

(71a) 

and 

^  =  Hi  cos(Otr  +  Hi  smoXr  . 

(71b) 

The  relations  between  Ei,  Hi  and  Eq,  Hq  are  explored  in  Appendix  C  with  the  result 
that 

_ f  ^  f 

El  =Eq  cosfer  +  jEo  sinfcr 

(72a) 

_ ft  _  f  _  // 

El  =—Eq  sin^r  +  ^o  coskr 

(72b) 

_  r  _ f  _ ff 

Eq  =  El  coskr  —  Ei  sin^r 

(72c) 

_  ff  _ f  _ ff 

Eq  -El  sinfcr+£i  cosfcr 

(72d) 

and  a  like  set  of  equations  relating  Hq  and  Hi. 

The  mathematical  manipulations  remain  unchanged  so  Equation  66a  could  just  as 
well  be  written  as 

_  _  /  _ /  _  / 

S-Nq  +Ni  cosloXj.  +  N2  sxnloXj. 

where 

(73a) 

Nq  =(Ei  xHi  +Ei  xHi  j/2 

(73b) 
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_ /  _  / 

with  similar  expressions  (from  Equations  66c  and  b)  for  Ni  and  N2  • 

The  integration  for  P  is  over  a  surface  and  does  not  involve  dr,  thus  sin  2©tr  and  cos 
2©tr  can  also  be  pulled  from  under  the  integrals  so  that  we  can  write 


/  /  / 

P  =  Mq  +Mi  cos2cotr  +  M2  sin2<i»r 

(74a) 

where 

Mq  =  f  Nqj.  r^sm6ddd(l> 

A 

(74b) 

with  equivalent  expressions  for  M,'  and  Mj'. 

Since  P  must  be  the  same  in  Equations  74a  and  68a,  they  may  be  equated  and  the 
following  relations  found: 

Mq  =  Mq 

(75a) 

r  ' 

Mi  =  Mi  cos2kr-M2  sm2kr 

(75b) 

/  ' 

M2  =  Ml  sm2kr  + M2  cos2^r 

(75c) 

Ml  =  MiCOs2kr  +  M2sin2kr 

(75d) 

M2  =-Mism2kr  + M2COs2kr  . 

(75e) 

As  Grimes  (References  14  and  15)  observed,  Equation  68a  [or  Equation  74a]  can  be 
written  in  the  form  of  Equation  42  or,  its  equivalent,  Equation  50,  or,  more  generally. 
Equation  54a  and  b.  One  way  to  do  this  is  to  set 


P  =  Mq  +  Ml  coslax  +  M2  sin 2ox  =  Afo[l  +  cos(2tot  -  5)]  +  M4  sin(26»  -  S)  (76) 


and  solve  for  M4  and  8. 
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From  Equation  76 

Ml  cos2<or  +  M2  sin 2ox  =  Mo[cos 5 coslax  +  sin S sin 2ax] 
+  M4[sin2n»cos5  -  cos2o»sin5] 
which  reduces  to  the  pair  of  equations, 


Ml  =  Mq  cosS-  M4  sin5 

(77a) 

M2  =  Mq sin5  +  M4 cos8  . 

(77b) 

Squaring  Equation  77a  and  b  and  adding  the  result  gives  (References  14  and  15), 

M4  =  ±(Mf  +  .  (78) 


The  simplest  way  to  find  5  is  to  divide  Equation  77a  by  77b  resulting  in 
Ml / M2  =  (Mq cos S - M4 sin S) /(Mq sin S  +  M^cosS)  . 

Factoring  out  cos5  from  the  numerator  and  denominator  on  the  right-hand  side  and 
solving  for  tan  5  gives 

tan5  =  (M0M2  - MiM4)/(MqMi  +  M2M4)  .  (79) 

Of  course,  if  retarded  time  had  been  used  instead,  one  would  have 

P  =  Afo[l  +  cos(26()f^  -  5')]  +  M4  sin(2<at;.  -  5')  (80) 


where 


Ma  =± 


M?  +M^  -mI 


n1/2 


V 


J 


(81) 
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and 


tanS'  =  (moM2  -  Ml  m/ +  M2  m/ j  . 


(82) 


Substituting  Equation  75d  and  e  into  Equation  81  gives 

M4  =  {±Mi  cos^  2kr + Mf  sin^  2kr  +  2M\M2  sin 2krcos2kr  +  cos^  2kr 

1/2 

+  Mi^  sin^  2kr  —  2M1M2  sin  2krcos2kr-MQ^  (83) 

=  ±(m?+m|-mJ)^^^  =  M4  . 


This  last  result,  taken  in  conjunction  with  Appendix  B,  shows  that  if 

(2  2 

Mf +  MI-M5I  is  chosen,  one  of  the  two  possible  solutions  of  this  general 


form,  as  discussed  in  the  previous  section  results,  whereas,  if  the  minus  root  is  chosen, 
the  other  solution  will  result.  On  the  other  hand,  if  one  works  with  retarded  time,  the 
same  two  solutions  would  result.  Using  the  plus  minus  subscript  to  indicate  which  one  of 
the  two  subscripts  in  Equation  83  is  chosen,  one  has 


/ 

2(iXr-S+  =2(Ot-5± 


which,  when  solved,  yields 


d±  =5±-2kr  .  (84) 


This  can  be  shown  by  direct  substitution. 

Thus  it  makes  no  difference  as  to  whether  retarded  time  is  used  or  not,  the  same  two 
possible  answers  of  this  same  general  form  result,  although  antenna  fields  are  generally 
most  easily  and  naturally  written  in  terms  of  retarded  time. 

If  Equation  76  is  to  be  interpreted  in  terms  of  resistive  and  reactive  power  flow,  only 
one  of  the  two  possible  forms  can  be  correct,  although  M,,  (for  a  one  port,  two  terminal 
antenna)  would  be  proportional  to  the  radiation  resistance  and,  without  a  priori 
knowledge  as  to  which  form  is  correct,  M4  could  at  least  be  taken  as  representing  the 
absolute  value  of  the  radiation  reactance. 
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Unfortunately,  this  interpretation  (Reference  22)  is  on  somewhat  shaky  physical 
grounds.  The  radiation  resistance  seems  to  represent  a  reasonable  physical  concept. 
Clearly  the  power  radiated  must  be  supplied  by  the  signal  source  driving  the  antenna. 
One  can  find  this  power  radiated  by  integrating  the  Poynting  vector  over  the  surface  of 
any  sphere  surrounding  the  antenna.  If  one  uses  the  complex  Poynting  vector,  one  takes 
(1/2)  the  real  part  to  find  the  time  average  power  radiated;  if  one  uses  the  time  domain 
form  of  the  Poynting  vector,  one  looks  for  that  part  that  is  independent  of  the  radius  of 
the  sphere  chosen.  If  the  power  radiated  is  given  in  terms  of  the  square  of  the  input 
current  to  the  antenna  terminals,  then  the  coefficient  of  that  square  is  the  radiation 
resistance,  and  if  there  are  no  other  losses,  such  as  conduction  losses  in  the  antenna,  this 
must  represent  the  resistive  part  of  the  input  impedance  as  seen  at  the  antenna  terminals. 

The  reactive  power  concept  is  not  so  clear,  M4,  the  reactive  power  flow  depends  on 
the  radius  of  the  sphere  chosen,  with  M4  generally  increasing  with  decreasing  r,  since  the 

terms  in  S  and  ^  involved  with  the  reactive  power  flow  go  as  1/r^,  1/r^,  etc.  The  sphere 
that  just  encloses  the  antenna  gives  the  largest  reactive  power  flow  and  is  closest  to  the 
antenna  terminals.  However,  as  seen  in  Figure  1,  unless  the  antenna  completely  fills  this 
minimum  sphere,  there  remains  considerable  volume  exterior  to  the  antenna  for  energy 
storage  and  thus  reactive  power  flow  within  the  sphere  but  exterior  to  the  antenna. 
Although  it  would  probably  have  to  be  done  numerically,  the  Poynting  vector  could  be 
integrated  over  a  surface  more  closely  conforming  to  the  antenna's  physical  shape. 
However,  as  discussed  in  Appendix  A,  the  theoretical  fields  close  to  an  antenna  are  of 
doubtful  correctness  in  terms  of  the  usually  given  electromagnetic  solutions,  and 
integrating  using  incorrect  fields  would  likely  add  more  error  than  correction. 

In  the  case  of  Figure  5,  a  capacitively  loaded  dipole,  the  capacitors  add  reactance  not 
included  in  the  contributions  of  the  external  fields  and  thus  there  is  always  some 
reactance  due  to  antenna  geometry  not  included  in  M4. 
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FIGURE  5.  A  Capacitively  Loaded  Dipole. 


Regardless  of  these  caveats  on  the  use  of  M4  to  find  radiation  reactance,  there  seems 
to  be  an  even  bigger  problem  in  recasting  Equation  68a  in  the  form  of  Equation  76  and 
interpreting  it  in  terms  of  resistive  and  reactive  power  flow.  To  compare  Equation  76 
with  the  circuit  analysis  of  the  previous  section  requires  that  M4  be  real,  which  in  turn 
requires  that  Mj^  +  ^  in  Equation  78,  otherwise  M4  is  imaginary.  In  Reference 

41  it  is  shown  that  for 


E  =  ±jnH  ,  (85) 

M,  and  Mj  are  zero.  Fields  meeting  this  condition  are  known  as  Rumsey  fields 
(References  41  and  42)  and  are  not  especially  esoteric.  Circularly  polarized  fields  have 
Ais  property. 

If  these  radiation  reactance  (and  resistance)  concepts  have  problems,  perhaps  one 
should  instead  deal  with  the  stored  energy  in  the  antenna  fields,  and  from  this  find  an 
antenna  "Q,"  and  an  operating  bandwidth.  This  approach  will  be  discussed  in  the  next 
section. 


ANTENNA  Q  AND  BANDWIDTH 


The  basic  definition  of  the  quality  factor,  Q,  is  given  by  (Reference  29) 
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Q  _  6)(energy  stored  in  the  circuit)  _  ;r(energy  stored  in  the  circuit) 
average  power  loss  energy  loss  per  half  cycle 


(86) 


Written  out  symbolically 


Q  =  G)UIWi  . 


(87) 


The  instantaneous  energy  stored  in  inductors  and  capacitors  is  found  from  the 
expressions 


ULi  =  U^/2 

(88a) 

Uci  =  CV^I2 

(88b) 

where  the  subscript  i  indicates  instantaneous  value  since  I  and  V  are  generally  functions 
of  time. 

The  power  dissipated  in  a  resistance  is 


WLi  =  IV  =  fR  =  V^/R  .  (89) 

Equation  86  specifies  the  "time  average  power  loss."  If  the  current  through  the  resistor  is 
of  the  physically  realizable  form 


I  =  Iq  cosax 


(90) 


then 


Wii  =  RIq  cos^  ox  =  RIo{l  +  cos2taf)/2 


(91) 


The  time  average  power  loss  is  then 

Wl  =  R1^/2.  (92) 

If  the  current  is  assumed  to  be  in  the  physically  non-realizable  time  harmonic  form 
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(93) 


and 


Wl  =  R1I*/2  =  RI^/2 


(94) 


the  time  average  power  loss  has  been  found.  One  could  recover  W^.  from  the  time 
harmonic  complex  notation  form  by  employing 


Wu  =  R{II *  +Re//)/2  =  (i?/2)Re(// *  +77)  . 


(95) 


Running  the  same  current  through  an  inductor  in  series  with  this  resistance,  the  stored 
energy  is 


Uu  =  {L/2)I^ cos^  ajt  =  7^1(1  +  cos2mf)/4  . 

Applying  this  directly  to  Equation  87  results  in  a  time  dependent  Q 


(96) 


Qt  = 


cbLj} +  cos2ax) 
2R 


(97) 


and  while  there  is  nothing  wrong  a  priori  with  this,  Q  is  usually  given  as  a  time 
independent  quantity.  Curiously  enough,  if  one  changes  the  definition  of  Q  in  Equation 
86  in  the  denominator  from  being  the  time  average  power  loss  to  instantaneous  power 
loss,  one  would  have 


Qn  = 


£b7o  7-(l  +  cos  2a)r)  /  4 
7?7o(l  +  cos2far)/2 


=  COU2R 


which  is  time  independent. 

If  U  is  defined  to  represent  the  peak  energy  stored,  which  is 


U  =  iIl/2 


(98) 


(99) 


one  finds  that 
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Q  =  ((h/2l/2)/(/?/J/2)  =  a£/R 


(100) 


This  is  the  usual  equation  given  for  the  Q  of  an  inductor  and  a  series  resistance,  and  it  is 
in  this  sense  that  many  engineering  texts  define  Q  (Reference  29). 

If  one  were  to  use  the  average  energy  stored  rather  than  the  peak  energy,  one  would 
have 


Via  =  lUl*  (101) 

and 

Qa^diLllR  ,  (102) 


which  is  the  same  as  in  Equation  98. 

It  is  this  definition  of  Q  that  is  used  by  other  texts  (Reference  43).  Again,  =  Q/2, 
but  this  is  not  always  the  case  and  this  will  be  discussed  in  detail  later.  In  this  case  we 
also  have  =  Q^,  but,  in  fact,  is  not  generally  independent  of  time  as  shall  be  seen. 

Suppose  the  series  inductor  is  replaced  with  a  series  capacitor.  The  time  harmonic 
voltage  across  the  capacitor  is  given  by 


y/  =  IcoC 


(103) 


so  that 


Vc  =  Rejy^.  j  =  [Iq! coC)smm  . 


(104) 


The  energy  stored  in  the  capacitor  is 

Uci  =  (c/q  /to^C^lsin^ cm  =  {Jq  /4aP'Cj{l  -  cos2cm)  . 


(105) 


The  peak  energy  stored  is 


Uc  =  lll2ap-C 


(106a) 
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and  the  time  average  energy  is 


Uca  =  Io/^0)^C  .  (106b) 


In  terms  of  peak  energy,  then,  the  Q  is 


e  =  (/o  /  2ct)2c)  /(riI  /  2)  =  1  /  G)CR  (107) 


which  is  the  usual  equation  for  the  Q  of  a  series  capacitance  and  resistance  (Reference 
29). 

In  terms  of  the  time  average  energy  stored,  the  Q  is 


Qa  =  l/2G>CR  .  (108) 


Thus  =  Q/2,  but  again  this  is  not  always  the  case. 

One  can  also  see  from  Equations  100  and  107  that  in  these  two  cases,  Q  =  1X1  /R,  i.e., 
the  absolute  value  of  the  reactance  divided  by  the  series  resistance.  Again,  this  is  not  true 
in  general. 

In  this  case,  it  is  important  to  note  that  is  given  by 


a,=[/o(i  -cos2tat)/46)C]/[R/o(l+cos2taf)/2]  =  tan^t»/2(oCR  .  (109) 


From  this,  it  appears  that  Qtt  might  not  be  a  useful  concept,  since  it  ranges  with  time 
from  zero  to  infinity. 

Turning  to  the  series  RLC  circuit  of  Figure  3,  one  could,  of  course,  directly  use  the 
previous  equations  by  using  the  current  as  given  in  Equation  90.  However,  Figure  3 

shows  an  applied  voltage  of  V^coscot,  so  it  is  instructive  to  use  Figure  3  directly  with  the 
given  applied  voltage,  rather  than  replace  it  with  the  current  of  Equation  90.  The  answer 
will  be  the  same. 

The  time  harmonic  current  is  given  by 
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r  =  V'IZ  =  Voej^  /[/?+ j{coL-l/ct)C)] 

=  exp{-jtan~^  X/ 
or 

r  =  (Vo  /|Z|)exp[7(Q»  -  <j))]  (110) 

where 

Z  =  and  0  =  tan“^  X/R  . 


The  time  dependent  current  is  then 

I  =  (Vo/\Z\)cos{qx-<I>)  .  (Ill) 

The  energy  stored  in  the  inductance  is  then 

Uu  =  (Lyo^/4|Z|2)[l  +cos2(6>r-^)]  (112) 

while  the  power  lost  in  the  resistance  is 

Wu=(viR/2\Z\'^)[l  +  cos2{cot-(l>)]  (113a) 

with  a  time  average  power  loss  of 

Hi,  =  V^fi/2|Z|2  ,  (U3b) 

The  time  harmonic  voltage  across  the  capacitor  is 

Vc  =  {JXcVq  /|Z|)  exp  j{(m  -  ^)  =  ( /  ojC|Z|)  exp[;(mt  -<p-7t/2)]  (114a) 

and  the  time  dependent  voltage  is 
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Vc  =  (Vo  /C()C1Z|)cos(g»  -  0  -  /r/2)  =  (VJ,  /  fi)ClZ|)sin(G»  -  0) 


(114b) 


The  energy  stored  in  the  capacitor  is  thus 


Uci  =  (cVo^  /26)^C^|Zp)[l  -  cos2((Br  -  ^)] 


(115) 


and  the  total  stored  energy  is 


Ui  =  Uu  +  Uci  =  (Vb^  /4lZp){l[l  +  cos2(ta  -  0)] + (l/0^c)[l  -  cos2(d»  -  0)]} 


Ui  =  |vp/4fi)lZp)[fiiL+l/taC+(taL-l/G)C)cos2(taf-0)] 


(116) 


The  time  average  stored  energy  is  then 


Ua  =  (v^  /4tt)lZp)(taL  +  1/toC) 


(117) 


The  peak  stored  energy,  however,  depends  on  whether  ©L  >  1/©C  or  not.  Thus 


/  -  o\f  taL  ;  £QL>l/n)Corn}>fi)o 
C/=  Vo/2<»Z  H 

\  /h/mC  ;  taL<l/n)Corn)<©o 


(118) 


where  the  resonant  frequency  is  given  by 


©0  ~  1  /  LC 


(119) 


Using  the  peak  energy  definition  of  "Q,"  then 


Q  =  C0LIR  ;  ©>©o 


(120a) 


Q  =  l/o>CR  ;  ©<©o  • 

On  the  other  hand,  the  average  Q  is  given  by 


(120b) 
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Qa=aiLI2R-v\llcoCR  . 


(121) 


At  very  high  frequencies  (O)  »  cOg)  where  1/2g)CR  is  negligible,  the  average  is 
0)L/2R  or  one  half  the  peak  energy  Q  as  given  in  Equation  120a. 

At  very  low  frequencies  (O)  «  cOg)  where  (ioL/2R  can  be  neglected,  goes  to 
l/2ci)CR,  also  one  half  the  peak  energy  Q. 

On  the  other  hand,  at  resonance  when  co  =  cOg,  the  peak  energy  Q  is 

j2d  =  L/RVIC=(l/R)VI7C  (122) 


while 


OaQ  =  L/2R^IIc  +  ^JLC/2CR  =  {1/R)-4lIc  =  Qo  .  (123) 


Thus,  in  a  series  resonant  circuit,  ranges  from  Q/2  to  Q.  It  can  be  shown  that  these 
limits  hold  for  any  circuit,  starting  with  the  premise  that  stored  energy  is  always  positive, 
(Uj  >  0),  and  that  it  consists  of  a  time  dependent  part  that  varies  as  2(0  and  a  time 
independent  part. 

It  follows  that  a  stored  energy  versus  time  plot  must  look  like  Figure  6,  which 
represents  a  sinusoidal  waveform  plus  a  dc  component  such  that  the  function  is  always 
positive.  From  this  figure,  if  =  0,  then  U  =  U/2  and  =  Q/2,  as  is  the  case  for  a 
series  resistor  and  capacitor. 

On  the  other  hand,  if  U-  is  small  (the  time  varying  part  is  small),  =  U  and  Q,  = 
Q.  At  resonance  in  the  series  resonant  circuit,  the  time  varying  part  of  the  stored  energy 
is  zero  (the  total  reactance  is  zero)  and  =  Q.  Thus  it  is  clear  that 


QI2<Qa^Q 


(124) 


for  all  circuits. 
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FIGURE  6.  A  Typical  Stored  Energy  as  a  Function  of  Time  Plot. 

So  far  the  examples  covered  have  been  series  circuits;  parallel  circuits  must  be 
discussed  also.  Consider  a  capacitor,  C,  in  parallel  with  a  resistance,  R,  across  which  is 

applied  a  voltage  VoCosmt.  The  instantaneous  energy  stored  is 

t'o=(cVo/4)(l  +  cos2£W) 
while  the  power  loss  is 

Wu  =  cos^ axIR^)  =  yo^(l  +  cos2ta)/2iR  . 

Then,  in  terms  of  peak  energy  stored,  the  Q  is 

Q  =  /2)/(Vo^  /2r)  =  OCR  (125) 


while  in  terms  of  time  average  energy  stored  it  is 


NAWCWDTP8419 


Thus  the  peak  energy  Q  is  the  inverse  of  the  series  RC  case,  but  the  time  average  Q^^is 
still  one  half  of  the  Q  as  calculated  from  the  peak  energy. 

Considering  the  case  of  an  inductance,  L,  in  parallel  with  a  resistance,  R,  the  current 
through  L  with  the  same  applied  voltage  is 

Il=(VQlo£)sm(Ot  (127) 

and  the  instantaneous  stored  energy  is 

=(L/2)(Vo^/n)^L2)sin2Q»  =  (7^/4ffl2L)(l-cos2u»)  .  (128) 

The  Q  is 

Q  =  RI(dL  (129) 

i.e.,  the  inverse  of  Equation  100  for  the  series  RL  case.  The  time  average  Q  becomes 

Qa  =  RI2(oL  .  (130) 

Consider  the  complete  parallel  RLC  circuit  of  Figure  4.  The  total  stored  energy  is 

Ui  =  Uu^Uci  =  VQ\ilaL+(oC+{(oC-\l(oL)co^2(oi\lAa>  .  (131) 

The  time  average  stored  energy  is 

Ua  =  Vo{o)C+\lcaL)IA(o  ,  (132) 

while  the  peak  stored  energy  is 


VoCn  ;  (aC^McoL 
|Vo^C/2®^L  ;  n)C<l/mL 

Thus 


(133) 
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Qa  =  R{(oC  +  \l(dL)l2 
and 


(134a) 


oaCR 

RIcoL 


;  coCtMoL 

;  Q)C<1/giL 


(134b) 


In  this  case,  after  comparison  with  Equations  120a,  120b,  and  121,  it  is  seen  that  the 
time  average  Q,  is  not  the  inverse  of  the  series  case,  rather  L  and  C  are  interchanged. 
Again,  at  resonance,  =  Q  and  is  in  accord,  as  it  must  be,  with  Equation  124. 

Reference  43  gives  two  equations  for  finding  the  Q  of  a  circuit.  They  are 


Qa  =  {(0l2R)3Kld(0  , 


(135a) 


Qa  =  {(0l2G)dBld(0  . 


(135b) 


We  have  added  the  subscript  a  to  Qa  because  it  will  now  be  shown  that  the  application  of 
either  of  these  equations  leads  to  the  time  average  Q^.  Consider  the  series  resonant  circuit 
of  Figure  3.  Now 


Qa  =  {col2R)d{(oL-\lG)C)ldco  =  {(oL  +  H  (0C)I2R 
which  is  of  course.  Equation  121. 

Now  apply  Equation  135b  to  the  parallel  circuit  of  Figure  4.  Then 

Qa  =  {aiRI2)d{(oC-\l(0L)ld(0  =  R(mC+l/mL)/2 

which  is  Equation  134a. 

Equations  135a  and  b  may  have  been  assumed  to  be  universal  but,  in  fact,  Equation 
135a  applies  only  to  series  circuits  and  135b  only  to  parallel  circuits. 

To  show  this,  consider  the  circuit  of  Figure  3.  The  input  impedance  is 

Z  =  R-\-j{aiL-\l(oC)  =  R+jX  (136a) 


and  the  input  admittance  is 
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r=G+;5  =  l/Z  =  {R-jX)/^R^  +  . 

Applying  Equation  135b, 


(136b) 


Qa  =  (fi)/2G)a[-X/(/?2  +X^)]/<9o 
or 

Qa  =  -  i?2  +  R^^dKIda 

and 

Qa  =  [(fi>L  +  l/Q)C)/2/j][(x2  -  +i?2)J  . 

This  is  not  the  same  as  Equation  121,  but  rather  is  modified  by  the  factor  (X^- 
R^)/(X^+R^).  For  ©  very  small,  this  factor  goes  to  unity  since  1/  ©C  »  R  and  1/  © 
C»  ©  L  and  this  equation  would  give  a  correct  value  for  Qj.  Also,  for  ©  very  large,  the 

same  is  tme.  Otherwise  this  expression  is  incorrect,  although,  curiously,  as  ©  -»  ©o,  X  -4 
0,  and  the  above  equation  reduces  to 

Qa=-{(OL  +  \I(OC)I2R 


i.e.,  the  negative  of  the  correct  value.  The  same  procedure  of  finding  the  input 
impedance  of  the  parallel  circuit  of  Figure  4  and  applying  Equation  135a  yields  similar 
results. 

Any  series  circuit  consisting  of  N,  resistors,  Nj  inductors,  and  Nj  capacitors  can  be 
replaced  with  Figure  3  where  R  is  the  sum  of  the  resistors,  L  the  sum  of  the  inductors, 
and  the  capacitance  C  can  be  found  from 


1/C  =  l/Ci+1/C2+ . +  l/C'iV3  • 


(137) 


A  parallel  circuit  containing  N,  resistors,  Nj  inductors,  and  Nj  capacitors  can  be 
represented  by  Figure  4  where  C  is  the  sum  of  the  capacitances,  and 


1/R  =  1/Ri+1/R2+ . +  l/%i 

with 


(138a) 
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1/L=l/Li  +  1/I^+ . +  1/L;V2  •  (138b) 


It  follows  that  Equation  135a  applies  strictly  only  to  series  circuits  and  Equation  135b 
to  parallel  circuits  and  neither  is  exact  for  general  series-parallel  circuits. 

It  is  common  to  give  the  Q  in  terms  of  a  resonant  frequency  divided  by  a  3-dB 
bandwidth,  i.e., 


Q  =  /o/A/o 


(139) 


In  fact  this  expression  is  used  so  much  that  it  is  easy  to  subconsciously  accept  it  as  a 
fundamental  definition  of  Q,  instead  of  in  terms  of  the  stored  energy  and  the  power  loss 
per  cycle.  Upon  reflection,  however,  it  is  immediately  clear  that  this  expression  is 
somewhat  meaningless  in  terms  of,  say,  a  simple  series  R  and  L  circuit  where  there  is  no 
resonance.  Thus  Equation  139  cannot  be  applied  to  every  circuit,  only  those  that 
resonate.  It  follows  that  if  an  antenna  is  operated  in  a  frequency  region  where  it  is  not 
self  resonant,  then  the  “Q”  of  this  antenna  can  be  treated  in  terms  of  Equation  87,  i.e., 
stored  energy  and  power  loss,  but  not  directly  in  terms  of  Equation  139.  That  is,  given 
the  Q  and  the  operating  frequency,  one  cannot  use  Equation  139  to  find  the  bandpass, 
since  the  antenna  is  not  resonant.  Antennas  that  are  small  with  respect  to  a  wavelength 
are  usually  in  this  situation  (References  44  through  47). 

If  it  is  assumed,  for  now,  that  resonance  means  a  real  input  impedance  (the  reactance 
goes  to  zero  at  f(,),  then  one  can  talk  about  the  “bandwidth”  only  in  terms  of  the  antenna 
plus  its  matching  circuit,  or  at  least  in  terms  of  the  antenna  plus  whatever  reactance  is 
added  externally  to  cancel  out  the  antenna  reactance.  For  example,  suppose  that  around 
the  operating  frequency,  f^,  our  particular  antenna  can  be  well  represented  by  a  frequency 
independent  resistance  in  series  with  a  frequency  independent  capacitance,  i.e.,  R  and  C 
of  Figure  3.  The  “Qp”  of  this  combination  is  given  at  resonance  by  Equation  107  with 

©  =  ©0-  Adding  a  series  inductor,  L,  so  that  L  =  1/  ©q^C,  the  input  impedance  at  f(,  is  just 
R.  It  should  be  noted  that  the  antenna  is  still  not  “matched”  in  the  classical  sense,  unless 
it  happens  that  R  =  Z^,  the  characteristic  impedance  of  the  generator  being  used  to  feed 
the  antenna. 

The  current  through  the  series  RLC  circuit  of  Figure  3  in  time  harmonic  form  is 


Ic  =  Vcl[R+ j{a)L  - 1  /  coC)]  ( 140) 

where  V  =  RefV^]. 

The  time  average  power  loss  in  the  resistor  is  given  by  Equation  94  and  is 


Wl  =  V^RH  \r^  +  (©L-  1/©C)2]  . 


(141) 


At  resonance  the  time  average  power  loss  is 
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Wlq  =  V^I2R 


(142) 


Since  the  time  average  power  loss  decreases  for  frequencies  on  either  side  of  co,,,  one 

can  look  for  the  3  dB  frequencies,  ©,  and  where  the  power  dissipated  is  down  by  a 
factor  of  2,  i.e.. 


»i,I,2  /Wlo  =  {vo  +(«>l,2i-  ]}(2/!/ v^)  =  1/2 


or 


=  1/2 


which  results  in 


(fi)l  2^  - 1  /  (0\2C'f  = 


(143) 


The  resistance,  R,  is  always  positive  so  that  in  taking  the  square  root  of  both  sides, 
coiL-McOiC^R  ;  0>2>03q  (144a) 


and 


\Ig)iC-coiL  =  R  ;  MoiCxoiL  otodiXOq 


(144b) 


Equations  144a  and  144b  can  be  solved  explicitly  for  and  cOj  giving 


0)2 


RC  +  ^Ir^C^+4Lc]/2LC 


(145a) 


and 


Oi 


rV+4LC 


I2LC  . 


(145b) 
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If 


ACD  =  ©2 


(146) 


then 


©0  /  A©  =  /o  /  A/  =  (OqL/R  =  Qp 


(147) 


where  Equation  147  is  also  Equation  139. 

Since  no  approximations  have  been  made,  this  expression  is  exact  for  a  series 
resonant  circuit. 

It  is  interesting  to  note  that 


A©2  =  ©2  -©0  =  C^>QRC+^|(OQR^C^  +4-2  /2 


(148a) 


while 


A©i  =  ©0  -  ©1  =  ©0  j^©o^C'  -  ^jcOQR^C^  +4  +  2 


/2  . 


(148b) 


Thus  the  resonance  curve  is  not  symmetrical  in  ©,  since  A©2  ^  A®,.  In  fact 


A©2  >  A©i 


(149) 


because 


■sj(olR^C^+4>2  . 


(150) 


For  high  Q  circuits  where  R  is  small  enough.  A®,  will  become  virtually  equal  to  A®2, 
however. 

Of  course  the  Q  involved  is  for  the  whole  series  RLC  circuit,  however,  the  antenna 
was  assumed  to  be  a  series  R  and  C  with  a  matching  coil  to  make  the  impedance  real  at 
the  operating  frequency/resonance.  Since  both  the  capacitor  and  inductor  store  energy, 
one  can  quickly  see  that  the  total  energy  stored  might  be  greater  than  for  the  antenna 
alone.  It  can  be  seen  that  adding  any  forther  reactive  elements  for,  say,  the  purpose  of 
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matching  the  input  impedance  to  a  characteristic  impedance  (or  any  other  purpose  for  that 

matter)  might  raise  the  Q,  and,  by  implication  narrow  the  bandwidth  Af  in  Equations  148 
and  140.  In  fact,  it  can  be  shown,  given  a  circuit  with  some  Q,  that  the  addition  of 
reactive  elements  in  either  series  or  parallel  carries  the  following  condition  on  the 
resultant  Q,; 


Qr^Q  .  (151) 


This  proof  is  carried  out  in  Appendix  D. 

In  the  example  just  covered,  the  Q,  of  the  matched  circuit  is  lAnCR  for  frequencies 
below  or  equal  to  the  resonant  frequency,  i.e.. 


Qr  =  Q  ;  /^/o 


(152a) 


and  the  addition  of  the  inductance  leaves  the  Q  unchanged.  For  frequencies  above 
resonance,  however, 


Qr  =  (oL!  R  =  G)l  (OqRC  =  (oQq!  (Oq  =  aP'Q!  (Oq 


(152b) 


That  is,  Q,  is  greater  than  the  Q  of  the  antenna  alone.  Of  course  it  can  be  argued  that  near 
the  operating  frequency,  f(,,  the  Q  doesn’t  change  much  for  a  high  Q  circuit  because  f/f^ 
must  be  near  unity  to  stay  within  the  3  dB  bandwidth. 

In  summary,  if  an  antenna  could  be  faithfully  represented  by  frequency  independent 
series  resistance  and  capacitance  (and  the  authors  know  of  no  such  antenna),  the  Q  as 
determined  from  the  radiated  power  and  stored  energy  could  be  used  to  determine  the 
operating  bandwidth  and  frequency  with  a  series  inductance  for  matching. 
Unfortunately,  this  has  not  been  shown  to  be  the  case  in  general.  For  example,  a  parallel 
inductance  could  be  chosen  as  a  matching  element  for  the  hypothetical  series  R  and  C 
antenna  in  Figure  7. 


FIGURE  7.  Using  a  Parallel  Inductance  as  a 
Matching  Element  for  a  Hypothetical  Series 
R  C  Antenna. 
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Before  dealing  with  this  case,  one  needs  to  derive  the  equivalent  of  Equation  148  for 
a  parallel  resonant  circuit.  While  it  seems  that  this  should  be  straightforward,  there  are 
conceptual  difficulties  to  overcome.  In  using  the  power  loss  of  Figure  3  at  resonance  and 
the  3  dB  frequencies  to  define  the  operating  bandwidth,  there  was  no  difficulty  because 
the  current  I  through  the  (radiation)  resistance  had  an  inverse  relation  with  the 
impedance,  Z.  However,  if  Figure  4  is  used  as  drawn,  maintaining  a  constant  voltage 

Vocoscot  across  the  circuit,  the  power  lost  (radiated)  will  be  independent  of  frequency  as 

Vq^cos^  ©t/R. 

Typically  in  portraying  experimental  work  the  signal  source  is  shown  as  a  voltage 
generator  in  series  with  a  series  impedance  Zg  so  that  actual  circuits  look  like  Figures  8a 
and  8b.  Figure  8a  can  be  treated  exactly  as  Figure  3  with  R  replaced  by  R,.  =  R  +  Z,,,  as 
long  as  Zq  is  real,  and  the  loaded  Ql  becomes; 


Ql  =  a£IRT  ;  f>fo  (153a) 


Ql  =  \/(0CRt  ;  f<fo  .  (153b) 


Figure  8b  is  different.  The  characteristic  impedance  Zg  absorbs  some  of  the  generator 
output,  but  as  Zg  goes  to  zero  the  power  lost  (radiated)  in  the  (radiation)  resistance,  R, 
becomes  independent  of  frequency.  One  could  replace  the  voltage  generator  with  a 
current  generator  having  a  characteristic  conductance  across  the  generator,  but  while  this 
simplifies  the  mathematics  of  the  situation  because  parallel  conductances  add,  the 
problem  remains.  The  two  representations  are  exactly  equivalent  as  is  shown  in 
Appendix  E. 


FIGURE  8(a).  Series  RLC  Circuit  With  Signal 
Source  Shown  as  a  Voltage  Generator  in  Series 
With  Impedance,  Zg. 
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Zo 


(b) 

HGURE  8  (b)  Parallel  RLC  Circuit  With  Signal 
Source  Shown  as  a  Voltage  Generator  in  Series 
With  Impedance,  Zg. 

Another  approach  is  to  examine  Figure  3  and  observe  that  iZI  is  a  Tninimnm  at 

resonance,  rising  to  infinity  as  f  -»  0  and  f  ^  The  admittance,  lYI,  then  has  a  typical 
resonant  shape  with  a  peak  at  f  =  fg. 

Looking  at  the  absolute  value  of  the  admittance, 

|F|  =  l/[i?^+(a!L-l/mCf]^^^  (154a) 

then 

Yq  =  1/^  (154b) 

and  at  the  3  dB  points 

|^,2|/|^|  = =1/a/2  (154c) 

which  results  in  the  equation 

{(Oi,2L-y(Oi2cf  =R^ 

which  is,  of  course,  just  Equation  143  and  the  solution  must  be  Equation  147  as 
previously  found. 
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In  Figure  4  it  is  the  absolute  value  of  the  impedance  that  is  maximum  at  resonance 
and  goes  to  zero  as  f  0  and  f  ->  «>.  The  admittance  of  Figure  4  is  just 


Y  =  G+j{oiC-\lo£)  (155) 

where  G  =  1/R.  Then 

Zo=l/G  (156) 

so  that  the  3  dB  relation  is 

|Zu|/|Zol  =  G/[g^  +(<i)i,2C-  ]'”  =  1/V2  (157) 

resulting  in  the  equation 

(<»l,2C-l/«)i,2if  =G^  .  (158) 


Equation  158  is  of  the  same  form  as  Equation  143  with  C  and  L  interchanged  and  R 
replaced  by  G.  Thus  the  solution  must  be  given  by  Equations  146  and  147  with  this 
interchange  and  substitution,  i.e. 


Affl  =  G/C  (159) 

and 

cOq!  Lo)  =  (Of,CIG  =  cOoCR  =  Qp  .  (160) 


A  parallel  resonant  circuit  has  the  same  relation  between  Q  and  the  resonant 
frequency  divided  by  the  3-dB  bandwidth  as  the  series  resonant  circuit. 

Of  course  not  all  circuits  are  wholly  series  or  wholly  parallel  but  may  be  some 
combination  of  series  and  parallel  elements. 

Figure  7  is  a  simple  example  of  such  a  circuit.  If  one  cares  to  analyze  this  circuit,  one 
will  find  it  far  more  algebraically  “messy”  than  a  simple  series  or  parallel  RLC  circuit.  If 
the  series  R  and  C  represent  some  “mythical  antenna,”  one  will  find  the  concepts  of 
matching  and  bandwidth  considerably  more  complicated. 
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CONCLUSION 


This  report  has  been  devoted  to  examining  some  of  the  basic  ideas  and  assumptions 
that  are  involved  in  the  theoretical  analysis  of  small  antennas.  Poynting’s  theorems  in 
both  the  time  and  frequency  domain  cases  have  been  derived  and  the  conditions  under 
which  they  are  mathematically  applicable  carefully  examined.  The  relationships  with 
stored  energy,  both  time  dependent  and  time  average,  were  derived. 

It  was  shown  that  the  concept  of  integrating  Poynting’s  vector  over  a  surface 
surrounding  an  antenna  to  find  the  power  flow,  although  valid,  is  not  generally  the  same 
as  or  even  a  part  of  Poynting’ s  theorem. 

Resistive  and  reactive  power  flow  into  circuits  was  examined  and  compared  with  the 
power  flow  across  a  closed  surface  surrounding  an  anteima.  It  was  pointed  out  that  while 
using  this  approach  to  find  the  input  or  radiation  resistance  of  a  single  port  antenna  can 
often  be  justified,  the  situation  as  regards  finding  antenna  input  reactance  is  on  very 
shaky  ground  indeed. 

Because  small  antennas  are  frequently  referred  to  in  terms  of  “Q”  and  bandwidth, 
these  concepts  were  explored  for  series  and  parallel  circuits.  Even  here  much  room  for 
misinterpretation  and  misconceptions  was  found.  Indeed  it  was  found  that  all  of  the 
references  quoted  here  dealing  with  limitations  on  the  performance  of  small  antennas 
(References  1  through  24)  could  be  criticized  on  the  grounds  of  misapplying  the 
fondamentals  covered  in  this  report. 

However,  much  more  remains  to  be  covered,  i.e.,  to  find  the  theoretical  limits  on 
antennas  that  are  small  with  respect  to  a  wavelength.  Some  of  this  work  has  actually 
been  done  in  the  form  of  notes  and  derivations  but  space,  time,  and  money  constraints 
prevent  its  inclusion  in  this  report.  Further  reports  are  planned  dealing  with  such  subjects 
as  some  particular  antenna  configurations  that  could  show  performance  much  improved 
over  the  Wheeler  limit  (Reference  2),  Rumsey  antennas  (Reference  42),  and  new 
applications  of  spherical  harmonics  and  multipoles  to  small  antenna  analysis. 
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Appendix  A. 

THE  “EXACT”  DIPOLE  ANTENNA  AND  ITS 
TANGENTIAL  BOUNDARY  CONDITION 


Consider  the  dipole  antenna  of  Figure  A-1,  constructed  of  conducting  wire  with  a 
diameter  of  a,  a  gap  width  g,  and  a  total  length  L.  If  the  wire  is  made  of  a  good 
conductor,  if  L  »  a,  the  gap  width  g  is  very  small  with  respect  to  a  wavelength,  and  L  is 
on  the  order  of  a  wavelength,  it  is  usually  assumed  (Reference  A-1)  that  Figure  A-2  is  a 
good  model  with  which  to  represent  Figure  A-1,  so  that  the  radiated  fields  can  be 
reasonably  found  and  calculated. 


L/2 


d 


cr 


L/2 


FIGURE  A-1.  Wire  Dipole  Antenna. 


The  model  in  Figure  A-2  could  be  said  to  be  benign  in  the  sense  that  increasing  the 
conductivity  of  the  wire,  decreasing  the  wire  diameter  a,  and  the  gap  width  g,  albeit  this 
limits  the  permissible  driving  voltage,  should  enable  the  real  antenna  to  arbitrarily 
approach  the  model. 

A  quite  general  solution  to  Maxwell’s  equations  can  be  found  using  a  vector  potential 
formulation.  For  the  time  harmonic  case,  this  vector  potential  is  given  by, 


57 


NAWCWDTP8419 


Alt  y>  K  -  ?■  1 


For  the  line  source  of  Figure  A-2,  this  reduces  to: 


L/2  f,-jkR 

A,=  J  m——dz' 


-L/2 


To  solve  the  dipole  problem  it  is  customary,  at  this  point,  to  assume  i 
distribution,  I(zO-  Typically  a  sinusoidal  current  distribution  of  the  form 


Hz')  =  /o  sin[it(L/2  -  |z'|)]  ;  -LI2<z'^LI2 

is  assumed. 


z 


FIGURE  A-2.  Model  for  Dipole  Antenna. 


(A-1) 


(A-2) 

current 

(A-3) 
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This  current  distribution  is  an  assumption  but  it  seems  a  reasonable  one,  with  the 
characteristics  of  a  standing  wave  and  going  to  zero  at  the  dipole  ends.  In  any  case,  the 
radiation  characteristics  of  a  dipole  do  not  appear  to  be  overly  sensitive  to  the  form  of  the 
assumed  current  distribution  (References  A-2  and  A-3),  with  assumed  rectangular, 
sinusoidal,  and  triangular  distributions  giving  very  similar  results. 

If  Equation  A-3  is  substituted  into  Equation  A-2  and  utilizes 


H  =  -VxA 


(A-4a) 


which  is  in  this  case 


1  dA^ 

H  dp 


(A-4b) 


one  obtains  terms  under  the  integral  that  are  exact  differentials  (Reference  A-1),  and 
carrying  out  the  integration,  the  magnetic  fifeld  is 


+e~^^  -2cos(ii:L/2)e"-'*''] 


(A-5) 


The  jionzero  electric  field  components  can  be  found  directly  from  this  by  taking  the 
curl  of  Jj  resulting  in. 


E  =  -  LI2)e~^^  IRi  +{z  +  LI2)e~j’^  IR2 

^  4;rp  t 

(A-6a) 


-2zcos(jtL/2)e”-^*''/r]  , 


IRi+e~^^^  IR2-2cos,{kLI2)e~^^''lr 


(A-6b) 


Within  the  limits  of  the  model  and  the  assumed  current  distribution,  these  appear  to 
be  exact  expressions  for  the  fields  everywhere  exterior  to  the  line  source.  It  is 
straightforward  to  transform  these  fields  from  cylindrical  to  spherical  coordinates  and 
make  the  assumption  r  »  L  and  find  the  usual  far  field  expressions  for  such  an  antenna, 
i.e.. 
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JIqb  cos{kLcosd)/2-cos{kL)l2 

Ha  = - 

^  27ir  sin0 


(A-7a) 


^0  ~ 


(A-7b) 


Since  no  dispute  of  these  results  has  been  found  in  the  literature,  it  is  assumed  that 
Equation  A-7a  and  b  represent  a  close  approximation  to  experimental  results  in  the  far 
field,  and  with  the  emphasis  on  near-field  measurements  in  Aese  past  several  years,  one 
might  expect  that  Equations  A-5  and  A-6a  and  b  have  been  experimentally  verified  to 
within  some  fairly  close  proximity  to  the  antenna  surface. 

Nonetheless,  Equation  A-6b  possesses  a  quality  which,  upon  discovery,  must  make 
any  engineer  who  started  his  career  in  metallic  waveguides  feel  at  least  uneasy.  Namely, 
the  tangential  electric  field  does  not  go  to  zero  at  the  antenna  surface,  or  even  anything 
close. 

This  is  most  easily  demonstrated  by  the  example  of  a  half-wave  dipole  where  L  =  X/2. 
It  then  follows  that  cos  (kL/2)  =  cos(2p/4>.)  =  cos  (7C/2)  =  0.  Next,  choose  an  observation 

point  P(p,z)  =  P(p,L/4)  =  P(p,l/8).  Choosing  z  =  0  would  simplify  the  math,  but,  for  p  = 
0,  this  corresponds  to  the  gap  on  the  real  antenna  and  enforcing  tangential  E  =  0  there 
might  be  considered  unreasonable. 

For  the  chosen  observation  point,  P,  Rj  and  Rj  become 


Rf  ={z-  L/2f  +  =  (A/8  -  A/4)^  +  =  A^  /64  +  (A-8a) 


=  (z  +  L/  2)^  +  p^  =  9A^  /  64  +  p^  (A-8b) 


=^^/64  +  p^  =R^  •  (A-8c) 


Equation  A-6b  for  E^  is  then 


=^^^{(A^/644-p2)  ^'^exp  -jk[X^ /64-\- 


(A-9) 


-f-|9A^/64  +  p^j  exp  -7k^9A^/64  +  p^j^ 


60 


NAWCWDTP8419 


If  one  now  takes  the  limit  as  p  goes  to  zero,  which  ought  to  be  tangential  E  for  a  very 
skinny  dipole,  there  results 


lim  =  2. 10834  Zq/q  /  TcX  •  (A-10) 

p-¥0 


It  could  reasonably  be  argued  that  even  though  tangential  E  must  go  very  close  to 

zero  for  the  real  antenna,  with  finite  conductivity  at  p  =  a/2  in  Figure  A-1,  it  is 
unreasonable  to  expect  this  for  the  line  source  in  Figure  A-2.  For  if  the  resistance  of  the 
wire  antenna  is 


Rr  =  a^lALa  (A-H) 


then  lim  R,  as  a  ^  0,  a  — >  is  indeterminate 

This  might  be  a  good  argument,  but  one  would  expect  that  tangential  E  ought  to  at 
least  have  a  local  minimum  somewhere  around  p  =  0.  This  doesn’t  seem  to  happen. 
Tangential  E  is  actually  maximum  at  p  =  0,  and  decreases  monotonically  after  that. 

To  show  this  one  can  take 


|£2p(;r^A^/4Zj/o)  =  F  =  l/(l  +  64/2)  +  l/(9  +  64/^) 


2  cos 

l9  +  64f-i 

^1  +  64/2^ 

/4) 

1  +  64/2^9  +  64/ 

(A-12) 


where  f  =  p/X,_ 

One  could  take  dF/df  and  show  that  it  is  always  less  than,  or  equal  to,  zero  and, 

therefore,  that  F,  and  hence  lEJ  decreases  monotonically  as  one  moves  along  p  away 
from  the  antenna.  It  is  easier,  however,  just  to  plot  F  versus  f  and  this  is  done  in  Figure 

A-3,  out  to  f  =  0.5  (a  half  wavelength).  F  is  maximum  at  p  =  0  and  decreases 
monotonically  with  f(p). 

Clearly,  tangential  E  does  not  only  not  go  to  zero  at  the  antenna  for  the  model  of 
Figure  A-2  and  the  assumed  sinusoidal  current  distribution,  it  doesn’t  even  go  to  a 
minimum. 
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Appendix  B 

ALTERNATE  FORMS  OF  THE  CIRCUIT  POWER  FLOW  EQUATION 
Suppose  we  have 

i(t)  =  Iqcos{qx-  P) 


The  power  flow  in  the  above  is  given  by 

Pit)  =  Vmt)  •  (B-1) 

Of  course  Vq  and  Iq  are  not  independent,  being  related  by  the  magnitude  of  the 

impedance.  The  phases  a  and  P  are  not  independent  either,  their  difference  being  a 
function  of  the  ratio  X/R  when  Z  =  R  +  jX. 

Carrying  out  the  operation  implied  in  B-1 

P(t)  =  IoVq  cos(QJt  -  P)cos{ctX  -  a) 
or 

P(t)  =  ^^^[cos(cf  -p)+ cos(2n»  -a-  p)]  (B-2) 

by  using  a  trigonometric  identity,  B-2  can  be  rewritten  as 

Pit)  =  [cos(a  -  /3)  +  cos(2tar  -ip  +  p- a)]  (B-3) 

or 
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P(t)  =  ^^^[cos(a  -P)  +  cos{P  -  a)cos(2<af  -  20)  -  sm{P  -  a)sin(2da  -  2^)] 

2 

and  finally 

P(t)  =  ■^^^[cos(a  -  P)[\  +  cos(2n»  -  2j8))  +  sin(a  -  j8)sin(2fi;f  -  2/5)]  .  (B-4) 

2 

Of  course  it  could  also  be  written  as 

P(0  =  •^^^[cos(a - P)  +  cos(2taf -2a  +  a- p)\  , 

2 

leading  to 

P(t)  =  ^^^^[cos(a  -  0){l  +  cos(2u»  -  2a))  -  sin(a  -  j8)sin(2u»  -  2a)]  .  (B-5) 

2 

As  seen  in  the  main  body  of  the  report  the  resistance  R  is  related  to  cos(a-P)  and  the 
reactance  X  to  the  sin  (a-P). 

The  exact  relation  can  be  found  by  observing  that 

7(0  =  Re^O^  ^  =  ^cos(q»  -  a  -  tan“^(X/P))  •  (®'6) 

R  +  jX  \Z\  ^ 

Thus 

/o=j^and/8  =  a  +  tan"^(X//?)  •  (2-7) 

Then 


R 


cos(a  -0)  =  cos|-tan  ^(X/P)]  =  cos^tan  ^(X/P)]  = 


(B-8a) 


and 
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sin(a-j8)  =  sin^-tan  ^(X/i?)j  = -sin|tan 


(B-8b) 


Both  B-4  and  B-5,  being  identities,  correctly  represent  the  power  flow,  however  only 
one  of  the  two,  as  explained  in  the  main  body  of  the  text  can  be  correctly  broken  up  so 

that  the  1+  cos  (2a)t  -  (j))  represents  resistive  power  flow  (where  (|)  is  either  2P  or  2a)  and 

sin  (2o)t  -  <{))  represents  reactive  power  flow.  Which  form  is  correct  requires  a  priori 
knowledge  about  the  circuit,  which  is  not  directly  found  in  the  integration  of  the  Poynting 

vector.  Nonetheless  the  multiplier  of  1  +  cos  (2mt  -  ()))  is  proportional  to  R,  so  R  can  be 
found  unambiguously  regardless  of  whether  the  proper  expression  is  chosen.  However, 

the  multiplier  of  sin  (2o)t  -  is  only  proportional  to  IXI,  so  that  only  IXI  could  be,  in  any 
case,  inferred  from  the  integration  over  the  time  domain  Poynting  vector.  In  any  case  B- 
4  and  B-5  are  of  the  general  form 


?  =  a[(1  +  cos(2£W  -  0))]  -1-  Bsin{2ox  -  (p)  . 


(B-9) 


Is  it  possible  to  write  either  of  these  (Equation  B-9)  in  the  form 


P  =  C[l  +  cos(2n»  -  0)]  -I-  Dsin{2QX  -  9) 


(B-10) 


where  C,  D,  or  0  are  different  than  A,  B,  or  (j)?  If  so  then  further  ambiguities  as  to  the 
value  of  R  and  X  (actually  IXI)  would  be  introduced  upon  the  integration  of  the  time 
domain  Poynting  vector. 

Set  B-10  equal  to  B-9  obtaining 


C[l  -I-  cos(2n)r  -  0)]  +  Dsin{2cot  -6)  =  A[1  +  cos(2mf  -  0)]  +  Psin(2<i»  -  0)  (B- 1 1 ) 


where  A,  B,  and  (|)  are  known,  C,  D,  and  6  are  unknown. 
Rewriting  B-1 1  as 


C  +  C  cos  0  cos  2G)r  +  Csin  0  sin  2for  +  D  cos  6  sin  2mf  -  D  sin  6  cos  2£0f 
=  A  +  A  cos  0  cos  2G)f  +  A  sin  ^  sin  2ox  +  Bcos  0  sin  2ox  -  P  sin  ^  cos  2(at 
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C+(Ccos0  -  Dsin0)cos2fi)f  +  (Csin0  +  Dcos0)sin  2cot 

(B-12) 

=  A +  {Acos(f>-  Bsin  0)  cos  2ctX  +  (A  sin  0  +  Bcos  (j))  sin  2(at 
To  be  true  for  all  values  of  t  one  must  have 

C  =  A  .  (B-13) 

This  can  easily  be  seen  by  setting  2ga  =  n,  then  2cot  =  -  Jt,  resulting  in  the  equations 

C  +  CcosO  -  DsinO  =  A  +  Acos<l>  -  Bsm(p 
and 

C-Ccosd  +  Dsin6  =  A-Acos<l)  +  Bsm(p  • 

Adding  these  equations  gives  C  =  A. 

Thus  B- 1 0  is  constrained  to  the  form 

?  =  A[l  +  cos(2taf-0)]+Dsin(2<a-0)  •  (B-14) 

Equation  B-12  becomes 

(A  cos  0  -  D  sin  0)  cos  2cjt  +  ( A  sin  0  +  D  cos  6)  sin  2ox 

=  ( A  cos  ^  -  B  sin  cos  2tor  +  ( A  sin  0  +  Bcos  sin  2ta 

To  be  true  for  all  values  of  t  one  must  have 

Acos0-Dsin0  =  Acos0-Bsin^  (B-15a) 

Asin0  +  Dcos0  =  Asin^  +  Bcos0  (B-15b) 

where  A,  B,  and  <j)  are  assumed  known  and  D  and  0  unknown. 
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To  eliminate  “D”  from  B-15a  and  B-15b,  multiply  through  by  cos 
resulting  in 

n 

Acos  6 - Dsindcosd  =  Acos0cos0-Bsin0cos6 
and 

Asin^0  +  Dsin0cos0  =  Asin0sin0  +  Bcos0sin0  . 

Adding  the  above,  one  obtains 

A  =  A(cos0cos0H-sin0sin0)  +  B(cos0sin0-sin0cos0) 
or 

A  =  Acos(0-^)  +  Bsin(0-0)  • 

Equation  B-16  can  be  solved  as  follows.  Rewrite  it  as 

A  ^  Btan(g-0) 

-\jl  +  taxP'{6  +  (l>)  'y/l  +  tan^(0-0) 

or 

A'^l  +  tan^(0-0)  =  A  +  Btan(0  -  ^)  • 

Squaring  both  sides  and  collecting  terms,  one  obtains 

(A^-B^)tan2(0-0)-2ABtan(0-0)  =  O  • 

Thus  either 

tan(0  -  0)  =  0 
or 

,,  2AB 


and  sin  6 


(B-16) 


(B-17) 


(B-18a) 


(B-18b) 


69 


NAWCWDTP8419 


and  the  solutions  to  B-18a  are 

d  =  <j)  »  (B-19a) 

d  =  (f>  +  n:  >  (B-19b) 

d  =  (f)-jc  •  (B-19c) 

Substitution  of  B-19a  into  B-15a  and  15b  leads  immediately  to  the  trivial  solution 

D  =  B  .  (B-20) 

Substitution  of  B- 19b  into  B-15a  and  15b  gives 

-Acos0  +  Dsin^  =  Acos0-Bsin0  (B-21a) 

and 

-Asin0-Z)cos^  =  Asin^-Scos^  .  (B-21b) 

Multiplying  the  first  equation  by  sin<|)  and  the  second  by  coscj)  and  subtracting  gives 

D  =  -B  .  (B-22) 

However  substitution  of  this  solution  back  into  B-21a  and  B-21b  results  in  A  =  -A  which 
is  not  allowed  unless  A  =  0  which  is  not  a  general  case. 

Substituting  B-19c  into  B-15a  and  B-15b  gives 

-A  cos  0  +  Dsin  0  =  A  cos  0  -  Bsin  0  (B-23a) 

-Asin0-2)cos0  =  Asin0  +  Bcos0  .  (B-23b) 

Again,  multiplying  the  top  equation  by  sin0,  the  bottom  by  cos0  and  subtracting  results  in 
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D  =  -B  , 


(B-22) 


but  with  the  same  result,  A  =  -A,  i.e.,  A  =  0.  Thus  B-19b  and  B-19c  are  not  general 
solutions. 

This  leaves  B-18b  as  the  only  other  possible  general  solution.  In  fact  since  it  is 
known  from  B-4  and  B-5  there  are  two  such  solutions,  it  is  the  other  general  solution. 
However,  to  find  D  in  terms  of  B  by  direct  substitution  of  B  18b  into  B-15a  and  B-15b  is 
algebraically  messy.  It  is  easier  to  return  to  B-15a  and  B-15b  and  rewrite  them  as 


A(cos0  -  cos^)  =  Dsin  0  -  Bsin  0 


(B-24a) 


and 


A(sin 6  -  sin 0)  =  Bcos0  -  DcosO 


(B-24b) 


Eliminating  A  by  multiplying  the  top  equation  by  (sin0  -  sin(j))  and  the  bottom  equation 
by  (cos0  -  COS0)  and  subtracting  gives 


Z)(l  — sin0sin0  — cos0cos0)  +  fi(l-sin0sin0  — COS0COS0)  =  0  • 


Either 


D  =  -B 


(B-25) 


or 


l-sin0sin0-cos0cos0  =  l-cos(0-0)  =  O  •  (B-26) 


Equation  B-26  is  the  (trivial)  solution  already  given  in  B-19a 

Therefore  B-25  corresponds  to  the  solution  given  in  B-18b.  This  can  be  verified  by 
substitution  of  B-25  and  B-18b  into  B-15a  and  B-15b. 

Thus,  it  follows  that  Equation  B-9,  corresponding  to  Equation  B-4,  can  be  rewritten 
as 


P  =  A 


,  _1  2AB  ^ 

1  +  cos  2(Ot  -  0  -  tan  — ^ - s- 


( 


-  Bsin|^2ta  -  0  -  tan' 


■1  lAB  \ 
A^-B^. 


(B-27) 
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which  must  then  correspond  to  Equation  B-5,  and  these  are  the  only  two  possibilities  for 
this  term. 

There  are,  however,  two  more  forms  of  the  power  flow  equation  that  might  tempt  one 
to  break  them  up  into  separate  resistive  and  reactive  power  flow  components  incorrectly. 
Consider  Equation  B-9.  Let 


(jy^  =  (j)  +  7r 


(B-28) 


Then 


?  =  a[i  +  cos(26»  -  01  +  ;r)]  +  Bsin(26>f  -  0i  +  ;r) 
or 

P  =  A[l-cos(2ft»-0i)]-Bsin(2<ar-0i)  .  (B-29) 

If  one  lets  02  =  0  -  (fe  and  0,  differ  by  2k  and  are,  in  effect,  the  same  angle.  The 
same  is  essentially  true  for  02=  7C  -  0. 

However,  it  follows  that  since  Equation  B-10  can  be  written  as 

P  =  A[l  +  cos(2o>f-0)]-Bsin(2{Hf-0)  , 
by  setting 

ei=e  +  n:  .  (B-30) 

it  can  also  be  written  as 

P  =  A[l-cos(2fi)t-0i)]  +  Bsin(2taf-0i)  •  (B-31) 

Thus,  from  the  preceding,  the  circuit  power  flow  equation  can  be  written  in  the 
separate  forms 

P  =  A[1  +  cos(2a»  -  0)] + Bsin(2fi»  -  0)  (B-32) 
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P  =  A[1  +  cos(2g»  -  0)]  -  Bsm{2ox  -  6) 

P  =  a[1  -  cos(2fflir  -  0i)]  -  Bsmilax  -<j>i) 
P  =  A[1  -  cos(2o»  -  ©j)]  +  Psin(2£i»  -  0i) 

where 


(B-33) 

(B-34) 

(B-35) 


(l>l  =  <!>  +  7C ,  6i  =6  +  k  ,  and  0  =  0  +  tan 


2AP 

A^-P^ 


(B-36) 


All  of  these  equations  are  correct  for  the  total  power  flow,  but  only  one  (either  B-9  or 
B-10)  can  generally  be  broken  into  separate  resistive  and  reactive  power  flows. 
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Appendix  C 

THE  RELATIONS  BETWEEN  THE  FIELD  COMPONENTS 
IN  STRAIGHT  AND  RETARDED  TIME 


Electromagnetic  fields  in  the  complex  time  harmonic  form  can  be  written  as 


E  = 


(C-1) 


From  tMs  the  time  domain  (sinusoidal)  physical  forms  can  be  found  by  taking  the  real 
part  of  E  (or  H  as  the  case  may  be). 

However,  antenna  fields  usually  come  in  the  form 


E  =  E^e-j^e^^  =  (f/  -  jE^y^’  (C-2) 


where 


tj.=t  —  rlc  •  (C-3) 

Thus  it  is  simpler  and  more  logical  to  write  the  antenna  fields,  both  time  harmonic  and 
time  domain,  in  terms  of  retarded  time. 

The  relations  between  and  are  simply  derived,  starting  with  Equation  C-1,  i.e., 
E  =  Eoej^e-j^’'ej^  =  V**'  =  EieJ^'  •  (C-4) 

Thus 

Eoe>^  =  Ei  (C-5) 

or 
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^£()  -jEq  j(cosA:r  +  7sinfcr)  =  £1  -jE^ 
and  from  this 

£1  =  £0  cosfcr  +  £o  sin^:r 

_ //  — .  /r  _  / 

£1  =£0  cosfcr-£o  sinfer  . 

Going  the  other  way 

I = V“'  = 

or 

^£1  -  jEx  ^{coskr- jsinkr)  =  Eq  -JEq 
and 

_  /  _ /  _ // 

Eq  =  El  coskr-Ei  sin  kr 

^  ff  —  //  _ / 

£0  =£1  coskr  +  Ei  sinkr  . 

A  like  set  of  equations  can  be  derived  for  the  magnetic  fields. 


(C-6) 

(C-7a) 

(C-7b) 

(C-8) 

(C-9) 

(C-lOa) 

(C-lOb) 
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Appendix  D 

EFFECT  ON  Q  OF  ADDING  PURELY  REACTIVE 
ELEMENTS  TO  ANY  CIRCUIT 


Suppose  there  is  a  box  (antenna)  containing  some  combination  of  resistances, 
capacitances,  and  inductances  in  any  series,  parallel,  or  series  parallel  arrangement.  This 
box  has  some  known  Q. 

Suppose  there  is  some  other  box  containing  only  lossless  capacitors  and  inductors  in 
any  arrangement  which  is  to  be  attached  to  the  first  box  for  some  purpose,  matching  at 
some  frequency. 

There  are  two  ways  to  attach  this  box,  i.e.. 


and 


The  Q  of  Za  was  originally  determined  by  applying  some  voltage  V  across  it  and 
driving  some  current  I  through  it.  Where  black  box  Xj  is  joined  to  Za  in  either  the  series 
or  parallel  sense,  V  and  I  might  change  if  one  drives  the  circuit  with  the  same  generator 
although  they  must  retain  the  same  relative  phase  and  amplitude  since 


V  =  ZaI  . 


(D-1) 
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The  power  flow  into  A  is  determined  from 


P^VI^Za?  =  V'/Za  . 


(D-2) 


After  the  addition  of  circuit  Xj,  in  either  series  or  parallel,  the  exact  same  V  across 
Za,  and  hence  the  same  I  can  be  maintained  by  changes  in  the  magnitude,  phase  and/or 
characteristic  impedance  or  admittance  of  the  driving  generator.  For  example  consider  a 
worst  case  scenario  where  circuit  Xj  approaches  a  short  and  is  placed  in  parallel  with 
circuit  Za- 

If  the  initial  measurement  is  carried  out  as 


then  the  same  V  and  I,  can  be  maintained  as  follows: 


It  is  postulated  that  this  is  done.  Now 


(D-3) 


where  Wl  is  the  time  average  power  loss  and  V  is  the  peak  stored  reactive  energy.  Wl  is 
unchanged  by  the  addition  of  Zj  if  V  and  I  are  unchanged  since  Zi  is  lossless.  The 
reactive  energy  stored  in  Za  is  also  unchanged  and  is  of  the  form 


VA(t)  =  Kai  +  Ka2  cos(2q»  -  (1>a) 


(D-4) 
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There  is  a  time  independent  part,  Kai,  which  must  be  positive,  and  a  time  dependent 
part,  Ka2  cos  (loot  -  (t)A).  can  also  be  written  as  a  positive  number  by  adding  or 
subtracting  n  from  (t)A.  That  is  if  <|)a  =  +  K,  then 

cos(2£0f - <l> a)  =  cos(2<af -<I>a~^)  =  ~ cos(2£(» -(I>a)  '  (D-5) 

It  will  be  assumed  that  this  operation  is  always  carried  out  so  that  Ka2  is  positive. 
Because  the  stored  energy  can  never  be  negative,  one  must  have 

The  Q  of  circuit  A  alone  is  thus 


_0)(KAi  +  Ky^2) 


(D-7) 


When  circuit  Zi  is  added,  whether  it  is  in  series  or  parallel,  and  the  driving  generator  is 
changed  to  maintain  the  same  V  and  I  across  and  through  circuit  A,  it  will  also  have 
stored  energy  of  the  form 


V2(r)  =  j  +  Ki2  COs(2£Of  ~(pi) 

where 


(D-8) 


Kii>ICi2 


(D-9) 


Again,  (j)!  is  chosen  by  adding  7i  so  that  K12  is  positive. 
The  total  time  dependent  stored  energy  is  thus 

U(t)  =  Ui(t)  +  UA(t)  =  Ki  +  K2  cos(2cot  -  (f) 


—  K-ax  +  jSTii  +  Ka2  cos(2c3r  —  0^4)  +  ^12  cos(2g)?  — 


(D-10) 
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By  breaking  D-10  into  cos  2(0t  and  sin  2(Dt,  one  finds  that 

Kji^2  +  ^12  ^  =  ^^2  cos  ^  > 

K/^2  +  ^12  sin  01  =  ir2  sin  0  • 

It  follows  that 

^A2sin0A  +  ^l2sin0i  . 

^A2  COS0^  +  ^12  COS  01 

Squaring  Equation  D-1  lb  and  D-1  Ic  and  adding  gives 

kI  =  KI2  +  KI2  +  2K/i2Ki2  cos(^A  -  #l)  • 

Thus  the  total  time  dependent  stored  energy  is 

U(t)  =  +  Kii  +  '\/^2^+-^^^  +  2^^42^^2^os(0^4-^  cos(2fi»  - 

The  peak  energy  stored  is 

U  =  KAi+Kii+^Kl2+Kf2+2KA2Kn  “S(^A  "  h) 

and  the  resultant  Q  is 

®f^Al  +'\/^A2  +  ^1^2  +2Ar^2^12COs(0^  -0l)  J 

j2  =  -i - 4  . 

Wl 

Clearly 


(D-1  la) 
(D-1  lb) 
(D-llc) 

(D-12) 

(D-13) 

(D-14) 

(D-15) 

(D-16) 


80 


NAWCWDTP8419 


or  the  total  Q  is  minimum  in  those  cases  where 


(b-17) 


(D-18) 


Then 


Q  =  <<K:M+Kn^K^2-Kn)  .  ^ (D-19a) 

or 

Q  =  ^  ,  p.,5b) 

L 


If  Ki2  =  Kn,  as  it  does  for  a  simple  capacitor  or  inductor.  Equation  D-19a  becomes 


Q  =  +  ^  (£>.20) 

Wl 


which  is  the  case  for  the  simple  series  RLC  circuit  below  resonance  as  discussed  in  the 
body  of  this  report. 

Note  that  if  Equation  D-19b  happens  to  be  the  case  and  Ka2  =  Kai,  then 


Hi  ““ 


since 


*:„  >  K-n  a  ■  ff>-22) 


It  has  been  shown  that  Q  ^  Qo  for  any  two  terminal  circuit  when  any  lossless  circuit  is 
added  to  it  externally.  One  could  add  further  circuits  (lossless)  in  any  fashion  such  as 
below. 
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It  is  readily  seen  that  by  treating  Xj  and  in  the  dotted  box  as  the  original  circuit, 
adding  X2  in  this  case  (in  series),  results  in  a  new  circuit  with  a  Q  equal  to  or  greater  than 
that  of  the  original  circuit. 

Thus  given  a  circuit  with  some  Qo,  adding  external  lossless  reactive  elements  to  it  in 
any  series,  series-parallel,  or  parallel  fashion  results  in  a  new  circuit  with  a  Q  greater  than 
or  equal  to  the  original  Qo- 
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Appendix  E 

EQUIVALENT  CURRENT  AND  VOLTAGE  GENERATORS 


Consider  a  voltage  generator  of  characteristic  impedance  Zq  feeding  an  impedance  Z, 
which  may  be  any  series,  parallel,  or  series-parallel  circuit. 


FIGURE  E-1.  Voltage  Generator. 


It  is  desired  to  replace  the  voltage  generator  with  a  current  generator  of  characteristic 
admittance  Yo,  Z  is  unchanged,  driving  the  same  impedance,  Z,  so  that  the  current 
through  (and  thus  the  voltage  drop  across 


FIGURE  E-2.  Current  Generator. 
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The  two  circuits  are  indistinguishable  to  the  right  of  plane  nun'  if 


r=i  . 


(E-1) 


Now  let 


jox 


I  =  R 


(E-2) 


Zf)+Z 


Let  both  Z  and  Zq  be  complex,  thus  allowing  for  a  complex  characteristic  impedance 
although  most  texts  treat  Zq  as  real.  Thus, 


/  =  Re 


Vpej^e-j^ 


^{R  +  Rof+(X  +  Xof 


(E-3) 


where 


0  =  tan 


X  +  Xq 
R  +  Rq 


(E-4) 


or 


Vq  cos(o!y  -  6) 

^{R+Bof+{X  +  Xof 


(E-5) 


For  the  current  generator  case 

Ig  +  I'  =  Is  •  (E-6) 

From  the  circuit  in  Figure  E-2,  it  is  obvious  that  the  voltage  across  Yq  and  Z  must  be  the 
same,  i.e., 

Ig,Zo=IoZ  .  (E-7) 
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Thus  Ir 


^  Z  ^ 

- 7-  +  1 


=  Igc  where  the  subscript  c  means  complex  notation.  Thus 


_ 

Z  +  Zq 


(E-8) 


therefore  /'  = 

Z  +  Zq 


or 


r  =  Re 


^0^0 

gjtan- 

|«  +  fio') 

[x+V] 

(E-9) 


where 


-iX  +  Xn 

0'  =  tan  ^ - iV  * 

R  +  Rq 


(E-10) 


Finally, 


lo 

r  - _ 

Zo 

cot 

^  ■ 

-e'  +  tan“^| 

Xo'/*6')) 

(x  +  V] 

l2  / 

1  +(^+:5^o 

f 

(E-11) 


Equating  Equations  E-1 1  and  E-5,  it  can  be  seen  that  I'  =  I  if 


Zb=Zo'=:^ 

^0 


(E-12) 


=  io\Zo\  =  Vo 


(E-13) 
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and 

^  =  — tan”^Xo  / Rq  =— tan  ^Xq/Rq 
thus  Figures  E-3  and  E-1  are  indistinguishable  as  far  as  Z  is  concerned. 


FIGURE  E-3.  The  Correct  Current  Generator  to  Deliver  the  Identical  Current 
to  the  Load  as  the  Voltage  Generator  cos  (Ot. 


Note  that  a  complex  characteristic  impedance  requires  a  phase  shift  between  the 
equivalent  voltage  and  current  generators. 
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